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14.  ABSTRACT 

Purpose 

The  goal  of  this  project  was  to  test  the  hypothesis  that  innate  immunity  protein  peptidoglycan  recognition 
protein  3  (Pglyrp3)  protects  from  inflammation,  obesity,  and  metabolic  syndrome  by  promoting  a  healthy  gut 
microbiota. 

Scope 

The  results  from  this  project  will  increase  our  understanding  of  genetic  factors  that  promote  obesity  and 
metabolic  dysfunction  and  the  role  of  gut  bacteria  in  the  development  of  metabolic  disease.  These  studies  will 
facilitate  future  development  of  novel  approaches  for  prevention  and  treatment  of  obesity  and  associated  diseases 
and  may  include  use  of  Pglyrp- controlled  microbiota  to  restore  and  maintain  beneficial  bacteria.  Such  novel 
approaches  are  essential  in  order  to  improve  the  health  of  military  personnel,  veterans,  and  their  family  members. 

Major  findings 

Pglyrp3'''  mice  on  high  fat  diet  (HFD)  gained  significantly  less  body  mass  and  had  significantly  less  adipose 
tissue  and  steatosis  than  WT  mice  on  HFD.  Pglyrp3 _/"  mice  on  HFD  had  lower  serum  triglycerides  and  developed 
significantly  less  glucose  intolerance  than  WT  mice  on  HFD.  We  further  determined  that  the  gut  microbiota  is 
altered  in  Pglyrp3~''  mice  on  HFD  compared  with  WT  mice  on  HFD. 

Up-to-date  report  of  results  and  significance 

Our  results  demonstrate  that  innate  immunity  proteins  modulate  the  development  of  obesity  and  metabolic 
disease  in  mice.  Pglyrp3  predisposes  to  diet-induced  weight  gain,  hyperlipidemia,  increased  adiposity,  steatosis, 
and  development  of  glucose  intolerance.  To  determine  the  role  of  the  gut  microbiota  in  the  development  of 
metabolic  disease,  we  colonized  Swiss  Webster  WT  germ  free  mice  with  Pglyrp3~'~  HFD  microbiome  or  WT  HFD 
microbiome.  In  these  initial  experiments  we  did  not  get  a  difference  in  weight  gain  or  in  metabolites  between  the 
two  groups  of  colonized  mice.  In  future  experiments  we  will  use  BALB/c  germ  free  mice  to  establish  the  role  of  the 
Pglyrp3~'~  HFD  microbiome  in  development  of  metabolic  dysfunction  and  obesity. 
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1.  INTRODUCTION 


Low-grade  inflammation  and  changes  in  the  gut  microbiota  are  associated  with  obesity,  however 
the  source  of  inflammation  and  the  specific  changes  in  the  gut  microflora  that  result  in  obesity  are 
poorly  understood  (1-6).  Innate  immunity  is  the  first  line  of  defense  and  there  is  increasing  evidence 
that  components  of  the  innate  immune  system  modulate  the  development  of  obesity  and  associated 
diseases  (7-9).  Peptidoglycan  recognition  proteins  (Pglyrps)  are  innate  immunity  proteins  that  regulate 
inflammation  in  response  to  bacterial  peptidoglycan  and  maintain  a  healthy  microbiome  (10-14). 
However,  the  role  of  Pglyrps  in  obesity  and  metabolic  syndrome  is  unknown.  The  goal  of  this  project 
was  to  test  the  hypothesis  that  innate  immunity  proteins,  protect  from  inflammation,  obesity,  and 
metabolic  syndrome  by  promoting  a  healthy  gut  microbiota. 

2.  KEYWORDS 

Obesity,  metabolic  diseases,  diabetes,  innate  immunity,  microbiome 

3.  ACCOMPLISHMENTS 
Major  goals 


Specific  Aim  1:  To  determine  the  role  of  Pglyrps  in  susceptibility  to 
obesity  and  metabolic  syndrome 

Months 

IU 

Major  Task  1:  Establish  obesity  model  in  /^/yrp-deficient  mice 

Subtask  1 :  Obesity  experiments  with  Pglyrp- deficient  and  WT  mice, 
weigh  mice,  calculate  food  consumed,  establish  protocols  for 
biochemical  assays 

1-9 

Drs.  Gupta  &  Park 

Subtask  2:  Collect  stool  samples,  blood  and  organ  samples,  assay  for 
cytokines,  lipids,  glucose,  insulin  sensitivity  test,  histopathology 

3-12 

Drs.  Gupta  &  Park 

Subtask  3:  Identify  bacteria  in  stool  samples,  analyze  biochemical 
data  and  identify  changes  in  bacterial  populations 

3-12 

Drs.  Gupta  &  Park 

Milestone(s)  Achieved:  Role  of  Pglyrps  in  the  development  of  obesity 
and  metabolic  syndrome  is  established 

Tasks  accomplished 

Dr.  Park  had  to  return  to  South  Korea  and  Dr.  Nunez  was  hired  for  this  project. 

Subtask  1  is  completed:  For  obesity  experiments  Pg/yrp-deficient  ( Pglyrp3~'~ )  and  WT  mice  were 
placed  on  high  fat  diet  (HFD)  or  on  low  fat  diet  (LFD)  for  30  weeks  and  weight  of  individual  mice  was 
measured  each  week.  The  amount  of  food  consumed  was  measured  for  each  cage  and  the  amount  of 
food  consumed  per  mouse  and  per  day  was  calculated. 

Our  results  demonstrate  that  Pglyrp3  predisposes  mice  to  diet-induced  weight  gain.  Pglyrp3'''  mice 
maintained  on  high  fat  diet  (HFD)  gained  significantly  less  weight  than  WT  mice  on  HFD  for  the  entire 
length  of  the  experiment  (Fig.  1).  There  was  no  difference  in  weight  gained  between  Pglyrp3~'~  mice  on 
HFD  or  low  fat  diet  (LFD),  whereas  there  was  a  significant  difference  in  weight  gained  between  WT 
mice  on  HFD  and  LFD.  There  was  no  difference  in  the  amount  of  food  consumed  between  Pglyrp3- 
deficient  (1.98  g  food  consumed/day/mouse)  and  WT  mice  (2  g  of  food  consumed/day/mouse). 
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Figure  1.  Pglyrp3'lm  mice  on  HFD  gained  significantly  less 
weight  than  WT  mice  on  HFD.  PglyrpJ1'  and  WT  mice  on  HFD  or 
LFD  were  weighed  each  week  for  30  weeks.  The  results  are  means 
±  SEM  of  18-21  mice/group.  Significance  of  difference  by  ttest  for 
PglyrpJ1'  HFD  versus  WT  HFD,  *P  <  0.05;  WT  HFD  versus  WT 
LFD,  *P<  0.05. 


Subtask  2  is  completed;  We  collected  blood,  liver,  and  adipose  samples  for  cytokines,  lipids,  and 
glucose  assays,  for  histopathology  and  for  insulin  sensitivity.  We  have  also  collected  stool  samples  for 
microbiome  |3-diversity  analysis  and  for  colonization  experiments. 


We  have  performed  assays  for  cholesterol,  triglycerides,  glucose,  and  insulin  and  our  results 
demonstrate  that  Pglyrp3  predisposes  mice  to  diet-dependent  hypertriglyceridemia  and  development  of 
glucose  intolerance.  Pglyrp3' mice  on  HFD  had  significantly  lower  serum  triglycerides  than  WT  mice 
on  HFD  and  Pglyrp3'''  mice  on  LFD  had  significantly  lower  serum  triglycerides  than  WT  mice  on  LFD 
(Fig.  2A).  There  was  no  difference  in  serum  cholesterol  between  the  two  genotypes  on  HFD  or  LFD. 


Both,  WT  and  Pglyrp3~'~  mice  on  HFD  showed  impaired  ability  to  restore  blood  glucose  in  response 
to  a  bolus  of  glucose,  however  WT  mice  developed  more  severe  glucose  intolerance  than  Pglyrp3~'' 
mice  (Fig.  2B).  These  data  demonstrate  that  Pglyrp3  predisposes  mice  to  development  of  diet- 
dependent  and  diet-independent  metabolic  dysfunction. 


A  B 


Figure  2.  Pglyrp3  predisposes  mice  to 
development  of  metabolic  dysfunction.  WT 

and  PglyrpJ1'  mice  on  HFD  or  LFD  for  30  weeks 
were  assayed  for  (A)  serum  triglycerides  and  (B) 
glucose  tolerance.  For  the  glucose  tolerance 
test,  mice  were  injected  intraperitoneally  with  2 
g  glucose/kg  body  mass  and  blood  glucose  was 
measured  at  indicated  times.  The  results  are 
means  ±  SEM  of  7-12  mice/group.  Significance 
of  difference  by  ffesf  for  Pglyrp3~~  HFD  versus 
WT  HFD,  *P  <  0.05,  **P  <  0.005;  for  PglyrpJ 
LFD  versus  WT  LFD  +P  <  0.05. 


WT  mice  on  HFD  had  higher  levels  of  serum  insulin  than  Pglyrp3~''  mice  on  HFD  (Fig.  3A),  which 
suggests  the  development  of  diet-dependent  insulin  resistance  in  WT  mice.  However,  there  was  no 
significant  difference  between  WT  and  Pglyrp3~'' mice  in  an  insulin  resistance  test  (Fig.  3B). 

We  next  measured  serum  levels  for  cytokines  TNF-a,  IL-6,  IFN-y,  and  IL-1  (3,  however  all  cytokines 
were  below  the  detection  limit  in  both  groups  of  mice.  These  results  suggest  that  the  increased  weight 
gain  in  WT  mice  on  HFD  compared  with  Pglyrp3 mice  on  HFD  was  not  associated  with  increased 
cytokines. 
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WT  HFD  -m-  Pglyrp3'/_ 


Figure  3.  Pglyrp3'''  mice  on  HFD  have  lower 
serum  insulin  than  WT  mice  on  HFD.  WT  and 

Pglyrp3~'~  mice  on  HFD  for  30  weeks  were  assayed 
for  (A)  serum  insulin  and  (B)  insulin  tolerance  test 
(ITT).  For  the  ITT,  mice  were  injected 
intraperitoneally  with  0.4  U/kg  body  weight  of 
Humulin-R  insulin  (HealthWarehouse)  in  sterile  PBS 
and  blood  glucose  measured  at  15,  30,  60,  90,  and 
120  min  after  injection.  The  results  for  (B)  are 
means  ±  SEM  of  7-12  mice/group.  Significance  of 
difference  by  ttest  for  Pglyrp3~'~  HFD  versus  WT 
HFD,  *P<  0.05. 


Subtask  3  is  completed:  We  have  analyzed  the  biochemical  data  and  our  results  demonstrate  that 
Pglyrp3  predisposes  mice  to  diet-dependent  hyperlipidemia  and  development  of  glucose  intolerance. 
We  isolated  DNA  from  stool  samples  of  Pg/yrp3-deficient  and  WT  mice,  and  assayed  for  bacterial 
diversity  using  pyrosequencing  of  the  variable  regions  of  bacterial  16S  ribosomal  RNA  (rRNA)  genes. 
We  compared  the  sequences  and  analyzed  the  changes  in  bacterial  ecology  using  Quantitative  Insights 
Into  Microbial  Ecology  (QIIME;  http://qiime.sourceforge.net),  Bayesian,  PyNAST,  and  UniFrac.  In  our 
preliminary  analysis,  the  microbiota  of  WT  and  Pglyrp3~''  mice  on  HFD  were  significantly  different  in 
their  p-diversity.  Once  the  analysis  is  completed  we  will  submit  the  complete  dataset  to  NCBI. 


Major  goals 


Specific  Aim  2:  To  determine  the  role  of  /^/yrp-controlled 
microbiota  in  susceptibility  to  obesity  and  metabolic  syndrome 

Major  Task  2:  Establish  obesity  model  in  germ-free  mice  colonized 
with  Pglyrp-deficient  or  WT  microbiota 

Subtask  1 :  Collect  stool  samples  from  Pglyrp3''~  and  WT  mice 

3-6 

Drs.  Gupta  &  Park 

Subtask  2:  Gavage  WT  germ- free  mice  with  Pglyrp3''~  and  WT  stools 
and  after  stable  colonization  of  gut  microflora  keep  mice  on  high  fat 
or  low  fat  diet 

6-12 

Drs.  Gupta  &  Park 

Subtask  3:  Obesity  experiments  with  stool  transplanted  mice,  weigh 
mice,  calculate  food  consumed 

9-12 

Drs.  Gupta  &  Park 

Subtask  4:  Collect  stool  samples,  blood  and  organ  samples,  assay  for 
cytokines,  lipids,  glucose,  insulin  sensitivity  test 

12-15 

Drs.  Gupta  &  Park 

Subtask  5:  Analyze  biochemical  data,  histopathology 

15-18 

Drs.  Gupta  &  Park 

Milestone(s)  Achieved:  Role  of  Pglyrp-contr oiled  microbiota  in 
susceptibility  to  obesity  and  metabolic  syndrome  is  not  demonstrated 

18 

Tasks  accomplished 


Dr.  Park  had  to  return  to  South  Korea  and  Dr.  Nunez  was  hired  for  this  project.  Dr.  Nunez  worked 
on  the  project  for  1 6  months,  but  took  up  another  position  after  that  because  of  the  lack  of  long-term 
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funding.  We  then  hired  Dr.  Gurses  for  this  project. 

Subtasks  1-5  are  completed:  We  collected  stool  samples  from  Pglyrp3'''  and  WT  mice  maintained 
on  HFD,  placed  them  immediately  in  reducing  buffer,  homogenized,  and  filtered  the  homogenate.  We 
used  the  stool  homogenates  to  colonize  Swiss  Webster  WT  germ  free  mice  with  Pglyrp3~'~  or  WT 
microbiota  and  maintained  all  mice  on  HFD.  The  weight  of  individual  mice  was  measured  each  week 
and  our  results  demonstrate  that  there  was  no  difference  in  weight  between  mice  colonized  with 
Pglyrp3~''  or  WT  microbiota.  There  was  also  no  difference  in  the  serum  triglycerides  and  glucose 
between  mice  colonized  with  Pglyrp3~''  or  WT  microbiota.  We  did  not  assay  for  cytokines  (part  of 
Subtask  4)  because  there  was  no  difference  in  weight  gain,  serum  triglycerides  or  serum  glucose 
between  mice  colonized  with  WT  or  Pglyrp3~'~  microbiota. 

The  lack  of  difference  in  weight  gain  and  metabolites  between  mice  colonized  with  Pglyrp3' or  WT 
microbiota  could  be  due  to  several  potential  reasons,  including  the  strain  of  germ  free  mice  used  for 
these  experiments.  We  used  Swiss  Webster  germ  free  mice,  which  are  outbred  mice  and  most 
frequently  used  for  similar  microbiome  colonization  experiments.  However,  these  mice  are  considerably 
different  from  BALB/c  mice,  which  we  have  used  for  other  experiments.  Recently,  BALB/c  germ  free 
mice  have  become  available  and  we  will  repeat  the  colonization  experiments  using  these  mice. 


In  a  separate  ongoing  project  in  my  laboratory,  we  determined  the  effect  of  another  innate  immunity 
protein  Nod2  on  diet-dependent  obesity  and  metabolic  disease.  We  demonstrate  that  Nod2~l~  mice  on 
HFD  compared  with  WT  mice  on  HFD,  gained  significantly  more  weight  with  increased  adiposity  (Fig.  4, 
ref.  15).  NodZ1'  mice  on  HFD  developed  hypercholesterolemia,  hypertriglyceridemia,  glucose 
intolerance,  and  insulin  resistance  (Fig.  5,  ref.  15). 


Figure  4.  Nod2 r/"  mice  on  HFD  gain 
significantly  more  weight  than  WT 
mice  on  HFD.  (A)  Body  weight.  (B) 
Representative  mice  shown  for  increase 
in  size  and  visceral  adipose.  The  results 
are  means  ±  SEM  of  10-12  mice/group. 
Significance  of  difference  by  ttest  for 
NodZ1'  HFD  versus  WT  HFD,  **P  < 
0.0005;  NodZ1'  HFD  versus  Nod?1'  LFD, 
mP  <  0.0005,  and  WT  HFD  versus  WT 
LFD,  AP  <  0.05,  and  AAP  <  0.005. 


Figure  5.  NodZ1  mice  on  HFD  develop  hyperlipidemia,  hyperglycemia,  glucose  intolerance  and  insulin 
resistance.  WT  and  NodZ1'  mice  were  assayed  for  (A)  cholesterol,  (B)  triglycerides,  (C)  glucose,  (D)  glucose 
tolerance,  (E)  insulin  and  (F)  insulin  resistance.  (D)  Mice  were  injected  intraperitoneally  with  2  g  glucose/kg  body 
mass  for  the  glucose  tolerance  test  and  (F)  with  0.4  U  insulin/kg  body  mass  for  the  insulin  resistance  test,  and 
blood  glucose  was  measured  at  indicated  times.  The  results  are  means  ±  SEM  of  6-10  mice/group.  Significance 
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of  difference  by  ttest  for  Nod2~'~  versus  WT  (HFD  or  LFD),  *P  <  0.05,  **P  <  0.005,  ****P  <  0.00005;  Nod2~'~  HFD 
versus  LFD,  mP  <  0.0005;  and  WT  HFD  versus  LFD,  AP  <  0.05,  AAP  <  0.005. 

We  identified  the  bacterial  diversity  in  WT  and  Nod2 ''  mice  maintained  on  HFD  or  LFD  for  30 
weeks  using  pyrosequencing.  The  microbiome  of  all  four  groups  of  mice  (WT  HFD,  Nod2'''  HFD,  WT 
LFD  and  Nod2'1'  LFD)  separated  from  each  other  in  the  Principal  Coordinate  Analysis  (PCoA),  with 
microbiota  from  each  group  showing  distinct  and  significant  separation  in  their  (3-diversity  (Fig.  6A,  ref. 
15,  16).  We  detected  many  statistically  significant  effects  of  both  the  Nod2  genotype  and  the  diet  on 
bacterial  (3-diversity;  with  many  significant  changes  in  the  abundance  of  various  bacterial  groups  in  the 
stools  at  all  taxonomic  levels  (class  diversity  shown  in  Fig.  6B,  ref.  15,  16). 


B 

WT  HFD 


WT  LFD 


Nod2'/-  LFD 


Nod2'/-  HFD 


■  Clostridia  ■  Candidatus  saccharibacteria 

■  Bacteroidia  ■  Deltaproteobacteria 

■  Gammaproteobacteria  ■  Flavobacteriia 

■  Bacilli  ■  Actinobacteria 

■  Epsilonproteobacteria  ■  Erysipelotrichia 


Alphaproteobacteria 
■  Cyanobacteria 
Spirochaetia 
Deferribacteres 
Betaproteobacteria 


Cytophagia 
■  Thaumarchaeota 
Sphingobacteriia 
Mollicutes 


Figure  6.  P-diversity  in  stool  microbiota  in  WT  HFD,  Nod2'''  HFD,  WT  LFD,  and  Nod2'1'  LFD  mice.  (A)  PCoA, 
each  sphere  with  confidence  ellipsoid  corresponds  to  a  stool  microbiome  of  one  mouse.  (B)  Class  abundance 
expressed  as  %  of  total  stool  microbiota  is  shown  as  average  for  all  bacterial  classes  (top  graphs)  and  for  17  less 
abundant  classes  expressed  as  %  of  total  low  abundant  classes  (excluding  two  most  abundant  classes,  Clostridia 
and  Bacteroida)  (bottom  graphs);  *,  classes  with  significantly  (at  P<0.05)  increased  abundance  in  WT  versus 
Nod2~'~  mice  (on  HFD  or  LFD);  #,  classes  with  significantly  (at  P< 0.05)  increased  abundance  in  HFD  versus  LFD 
groups  (for  WT  or  Nod2~'~  mice). 


We  next  considered  whether  the  altered  gut  microbiota  in  Nod2'''  mice  on  HFD  is  sufficient  for  the 
development  of  diet  dependent  obesity.  We  colonized  germ-free  pregnant  WT  Swiss-Webster  female 
mice  with  microbiota  from  either  Nod2''~  HFD  mice  or  WT  HFD  mice  and  placed  mothers  and  pups  on 
HFD  (Fig.  7A,  ref.  15).  Pups  born  to  the  mothers  colonized  with  Nod2'l~  HFD  microbiota  were  more 
susceptible  to  HFD-dependent  weight  gain,  hyperlipidemia,  and  hyperglycemia  than  pups  colonized 
with  WT  HFD  microbiota  (Fig.  7B,  ref.  15).  In  conclusion,  microbiota  from  Nod2'‘~  HFD  mice  can  transfer 
to  WT  outbred  germ-free  mice  enhanced  susceptibility  to  HFD-dependent  weight  gain  and  metabolic 
dysfunction. 


These  data  demonstrate  that  innate  immunity  molecules  regulate  development  of  diet-dependent 
obesity  and  metabolic  disease.  These  data  further  demonstrate  that  the  host  innate  immune  system 
together  with  the  diet  regulate  the  gut  microbiota  and  changes  in  the  innate  immune  system  together 
with  a  high  fat  diet  induces  dysbiosis,  which  promote  or  protect  from  obesity  and  metabolic  disease. 
Additionally,  these  data  further  demonstrate  that  our  approach  of  using  colonized  germ  free  mice  is  a 
relevant  approach,  with  the  potential  to  identify  the  mechanism  by  which  innate  immune  system,  gut 
microbiota,  and  diet  regulate  the  development  of  disease. 
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Figure  7.  Intestinal  microbiota  from  NodZ1'  HFD  mice  increases  susceptibility  of  germ-free  mice  to  HFD- 
induced  obesity  and  metabolic  dysfunction.  (A)  Germ-free  pregnant  mice  and  pups  were  colonized  with  WT  or 
NodZ1'  HFD  microbiota  and  maintained  on  HFD.  (B-E)  Pups  were  monitored  for  (B)  weight  gain  and  after  8  weeks 
assayed  for  (C)  liver  cholesterol,  (D)  liver  triglycerides,  and  (E)  serum  glucose.  (B)  Means  ±  SEM  and  (C-E) 
individual  data  of  19-21  mice/group.  Significance  of  difference  by  ttest ,  *P<0.05  and  **P<  0.001,  NodZ  '  HFD 
microbiome  vs.  WT  HFD  microbiome. 


In  summary,  our  results  demonstrate  that  innate  immunity  molecules  modulate  development  of 
diet-dependent  obesity  and  metabolic  dysfunction.  Specific  innate  immunity  molecules,  including 
Pglyrp3  and  Nod2,  promote  or  protect  from  diet-dependent  obesity,  hyperlipidemia  and  hyperglycemia. 
Our  results  further  demonstrate  that  in  the  absence  of  Nod2  and  in  the  presence  of  a  high  fat  diet  the 
gut  microbiota  is  altered,  which  contributes  to  the  development  of  disease  in  NodZ1'  mice  (15,  16). 

Opportunities  for  training  and  professional  development 

The  current  project  has  provided  all  participants  (Gupta,  Nunez,  Gurses,  Dziarski,  Caluag,  and 
Kirby)  with  several  opportunities  for  professional  development  and  has  resulted  in  increased  knowledge 
and  skill  in  the  areas  of  metabolism  and  obesity  research.  It  has  also  provided  opportunities  to  train  two 
undergraduate  students  (Caluag  and  Kirby)  in  the  process  of  research  and  in  the  development  of 
technical  skills.  Both  students  have  graduated  and  are  continuing  to  work  in  biomedical  research. 

How  were  results  disseminated? 

We  have  demonstrated  that  the  innate  immune  molecule,  Pglyrp3,  promotes  the  development  of 
diet-dependent  obesity,  hyperlipidemia,  and  hyperglycemia.  We  have  further  demonstrated  that  in  the 
absence  of  Pglyrp3  and  HFD  the  gut  microbiota  is  altered.  We  will  continue  to  work  on  this  project  to 
demonstrate  the  role  of  the  Pglyrp3'''  HFD-microbiota  in  protecting  from  diet-dependent  obesity  and 
metabolic  disease.  We  will  use  our  data  as  preliminary  results  for  extramural  funding  and  we  will 
publish  our  complete  findings  in  a  peer-reviewed  journal. 

We  have  recently  published  our  results  on  the  role  of  the  innate  immunity  molecule,  Nod2  and 
Nod2-regulated  microbiota  in  modulating  development  of  obesity  and  metabolic  disease  ( Scientific 
Reports,  2017,  7:548,  doi:  1 0.1 038/s41 598-01 7-00484-2,  ref.  15  and  16,  and  submitted  in  the 
Appendix). 


Plans  for  the  next  reporting  period 

There  is  no  next  reporting  period. 
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4.  IMPACT 

Impact  on  the  development  of  the  principal  discipline(s) 

Our  results  identify  a  novel  role  for  the  innate  immunity  gene  Pglyrp3  in  the  development  of  diet- 
induced  obesity  and  metabolic  dysfunction  and  further  identify  a  novel  susceptibility  gene  in  the 
development  of  theses  diseases.  The  genetics  of  obesity  in  humans  is  poorly  understood  and  it  is  likely 
that,  similar  to  other  complex  diseases,  a  vast  number  of  genes  contribute  to  susceptibly  to  obesity. 

In  a  separate  ongoing  project  in  my  laboratory,  we  discovered  that  the  innate  immunity  protein 
Nod2  protects  from  diet-dependent  obesity  and  metabolic  disease.  We  demonstrate  that  Nod2'''  mice 
on  HFD  compared  with  WT  mice  on  HFD,  gained  significantly  more  weight  with  increased  adiposity, 
hypercholesterolemia,  hypertriglyceridemia,  glucose  intolerance,  and  insulin  resistance  (Figs.  4  and  5, 
ref.  15).  We  further  demonstrate  that  the  Nod2'1'  HFD  microbiome  is  altered  and  contributes  to  the 
increased  weight  gain  and  metabolic  dysfunction  (Figs.  6  and  7,  ref.  15,  16).  These  results  confirm  our 
hypothesis  of  a  role  for  innate  immunity  and  gut  microbiome  in  regulating  metabolism  and  also  validate 
our  experimental  approach  to  determine  the  molecular  mechanism  of  regulation. 

Our  research  also  has  the  potential  to  identify  bacteria  species  that  promote  or  protect  from  obesity 
and  metabolic  dysfunction.  In  the  future  such  bacteria  may  be  used  as  probiotics  for  the  prevention  of 
metabolic  diseases. 

Impact  on  the  development  of  the  other  discipline(s) 

Pglyrps  have  an  established  role  in  innate  immunity,  bactericidal  activity,  and  inflammatory 
diseases.  Our  new  data  suggest  that  in  addition  Pglyrps  modulate  metabolism  and  also  supports  the 
emerging  concept  of  cross  talk  between  immune  and  metabolic  pathways. 

Impact  on  technology  transfer 

Nothing  to  report 

Impact  on  society  beyond  science  and  technology 

Nothing  to  report 

5.  CHANGES/PROBLEMS 

Changes  in  approach  and  reasons  for  change 

Nothing  to  report 

Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 

We  performed  the  experiments  as  described  in  Specific  Aim  2,  Subtasks  3  and  4.  We  colonized  WT 
germ  free  Swiss  Webster  mice  with  microbiome  from  WT  or  PglyrpJ1'  mice  maintained  on  HFD, 
maintained  all  colonized  mice  on  HFD  and  assayed  for  weight  gain  and  metabolites.  However,  there 
was  no  difference  between  the  two  groups  of  mice  for  weight  gain,  lipids,  or  glucose.  The  lack  of 
difference  in  weight  gain  and  metabolites  between  mice  colonized  with  Pglyrp3~'~  or  WT  microbiota 
could  be  due  to  several  possible  reasons,  including  the  strain  of  germ  free  mice  used  for  these 
experiments.  We  used  Swiss  Webster  germ  free  mice,  which  are  outbred  mice  and  most  frequently 
used  for  similar  microbiome  colonization  experiments.  However,  these  mice  are  considerably  different 
from  BALB/c  mice,  which  we  have  used  for  other  experiments.  Recently,  BALB/c  germ  free  mice  have 
become  available  and  we  will  repeat  the  colonization  experiments  using  these  mice. 

The  original  plan  was  to  hire  Dr.  Park,  a  postdoctoral  fellow  in  the  laboratory,  for  the  obesity  project 
as  Dr.  Park  has  extensive  experience  in  research  on  animal  models  of  inflammatory  diseases. 

However,  shortly  before  the  start  of  funding,  Dr.  Park  returned  to  South  Korea.  We  then  hired  Dr. 
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Nunez  and  he  started  work  on  the  project  in  September  2015.  Due  to  the  lack  of  long-term  funds  Dr. 
Nunez  left  our  lab  in  December  2016.  We  then  hired  Dr.  Gurses  who  started  work  on  the  project  in 
June  17.  These  changes  in  personnel  significantly  impacted  the  time  line  of  the  project. 

Changes  that  had  a  significant  impact  on  expenditures 

Nothing  to  report 

Significant  changes  in  use  and  care  of  human  subjects,  vertebrate  animals,  biohazards,  and/or 
select  agents 

Nothing  to  report 

6.  PRODUCTS 

Research  Material 

Our  research  has  led  to  the  discovery  of  unique  animal  models  to  investigate  the  role  of  specific 
innate  immunity  molecules  and  the  gut  microbiome  in  the  development  of  diet-dependent  obesity  and 
metabolic  diseases. 

In  the  U.S.,  more  than  68%  of  all  adults  are  overweight  or  obese  and  more  than  70%  of  Veterans 
who  receive  VA  care  are  either  overweight  or  obese.  More  than  half  of  obese  individuals  develop 
obesity-associated  diseases  including,  type  2  diabetes,  and  cardiovascular  disorders  (17).  Though  the 
primary  cause  of  obesity  is  an  excess  of  caloric  intake  over  expenditure,  host  genetics  and  the 
environment  contribute  significantly  to  the  development  of  disease.  There  is  increasing  evidence  for  the 
role  of  host  immune  system  and  the  gut  microbiome  in  the  development  of  metabolic  disease  (1-9). 
However,  the  molecular  and  biochemical  basis  of  this  influence  is  poorly  understood.  Furthermore, 
specific  bacteria  that  promote  or  protect  from  diseases  remain  mostly  unidentified. 

Research  using  animal  models  such  as  our  Pglyrp3 and  Nod2~'~  mice  is  crucial  to  obtain  a 
comprehensive  understanding  of  the  causes  of  metabolic  diseases.  A  detailed  understanding  of  the 
molecular  mechanisms  of  obesity  and  metabolic  dysfunction  will  enable  future  development  of  novel 
approaches  to  their  prevention  and  treatment  in  our  veterans  and  our  general  population.  These 
approaches  will  likely  include  re-balancing  the  gut  microbiome  using  probiotics. 

7.  PARTICIPANTS  AND  OTHER  COLLABORATING  ORGANIZATIONS 
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Dipika  Gupta 

Project  Role 

PI 

Months  Worked  on 
Project 

18 

Contribution  to 

Project 

Designed  the  experiments,  performed  some  of  the  experiments,  analyzed  the  data, 
and  wrote  the  grant  application,  annual  report  and  final  report 

Funding  Support 

Indiana  University 

Name 

Roman  Dziarski 

Project  Role 

pi 

Months  Worked  on 
Project 

12 

Contribution  to 

Project 

Designed  the  experiments,  performed  some  of  the  experiments,  analyzed  the  data 

Funding  Support 

Indiana  University 
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Name 

Ivan  Nunez 

Project  Role 

Postdoctoral  Fellow 

Months  Worked  on 
Project 

16 

Contribution  to 

Project 

Designed  the  experiments,  performed  most  of  the  experiments,  analyzed  the  data 

Funding  Support 

Name 

Serdar  Gurses 

Project  Role 

Postdoctoral  Fellow 

Months  Worked  on 
Project 

6 

Contribution  to 

Project 

Designed  the  experiments,  performed  some  of  the  experiments,  analyzed  the  data 

Funding  Support 

Name 

Tiffany  Caluag 

Project  Role 

Undergraduate  student 

Months  Worked  on 
Project 

9 

Contribution  to 
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Performed  some  of  the  experiments 

Funding  Support 

Indiana  University 

Name 

Kori  Kirby 
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Undergraduate  student 
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9 
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Funding  Support 

Indiana  University 
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Genetics  plays  a  central  role  in  susceptibility  to  obesity  and  metabolic  diseases.  BALB/c  mice  are  known 
to  be  resistant  to  high  fat  diet  (HFD)-induced  obesity,  however  the  genetic  cause  remains  unknown. 

We  report  that  deletion  of  the  innate  immunity  antibacterial  gene  Nod2  abolishes  this  resistance,  as 
Nod2~'~  BALB/c  mice  developed  HFD-dependent  obesity  and  hallmark  features  of  metabolic  syndrome. 
Nod2~'~  HFD  mice  developed  hyperlipidemia,  hyperglycemia,  glucose  intolerance,  increased  adiposity, 
and  steatosis,  with  large  lipid  droplets  in  their  hepatocytes.  These  changes  were  accompanied  by 
increased  expression  of  immune  genes  in  adipose  tissue  and  differential  expression  of  genes  for  lipid 
metabolism,  signaling,  stress,  transport,  cell  cycle,  and  development  in  both  adipose  tissue  and  liver. 
Nod2~'~  HFD  mice  exhibited  changes  in  the  composition  of  the  gut  microbiota  and  long-term  treatment 
with  antibiotics  abolished  diet-dependent  weight  gain  in  Nod2~'~  mice,  but  not  in  wild  type  mice. 
Furthermore,  microbiota  from  Nod2~'~  HFD  mice  transferred  sensitivity  to  weight  gain,  steatosis, 
and  hyperglycemia  to  wild  type  germ  free  mice.  In  summary,  we  have  identified  a  novel  role  for  Nod2 
in  obesity  and  demonstrate  that  Nod2  and  Nod2- regulated  microbiota  protect  BALB/c  mice  from  diet- 
induced  obesity  and  metabolic  dysfunction. 


The  increasing  prevalence  of  obesity  is  a  major  health  concern  worldwide1.  Obesity  is  a  risk  factor  for  other  dis¬ 
eases  and  more  than  half  of  all  obese  individuals  develop  metabolic  abnormalities,  which  include  hyperglycemia, 
hyperlipidemia,  hepatic  steatosis,  and/or  hypertension.  Though  the  primary  cause  of  obesity  is  an  excess  of  caloric 
intake  over  expenditure,  genetics  and  the  environment  play  a  central  role  in  regulating  metabolism  and  in  the 
development  of  these  diseases.  In  addition,  the  gut  microbiota  contributes  to  energy  homeostasis  in  the  host  and 
changes  in  the  composition  of  the  gut  microbiota  can  contribute  to  the  development  of  obesity  and  metabolic 
disease2, 3.  In  turn,  host  genetics  and  environmental  factors  shape  the  microbiota  and  changes  in  either  factor  can 
alter  intestinal  bacterial  populations.  However,  the  genetic  and  environmental  regulation  of  the  microbiota  is 
poorly  understood  and  the  specific  changes  in  the  gut  bacteria  that  contribute  to  the  development  of  obesity  and 
metabolic  disease  are  not  identified. 

Mouse  genotype,  similar  to  humans,  greatly  influences  metabolism  and  different  strains  of  mice  differ  in  their 
response  to  high  fat  diet  and  sensitivity  to  metabolic  diseases.  BALB/c  mice  are  naturally  resistant  to  the  devel¬ 
opment  of  high  fat  diet  (HFD)-induced  obesity,  however  the  mechanism  behind  this  resistance  is  not  known4. 
There  is  increasing  evidence  for  the  role  of  innate  immunity  in  the  development  of  metabolic  disease5, 6  and  in 
this  study  we  investigated  the  role  of  the  pattern  recognition  receptor  Nod2  in  the  resistance  of  BALB/c  mice  to 
HFD-induced  obesity 

Nod2  is  involved  in  innate  immunity  to  bacteria  and  is  expressed  in  pro -inflammatory  and  antigen  presenting 
cells,  such  as  monocytes  and  dendritic  cells,  and  also  in  structural  cells,  such  as  hepatocytes,  preadipocytes,  and 
intestinal  epithelial  cells7-10.  Nod2  is  a  cytoplasmic  sensor  for  bacterial  peptidoglycan  fragments,  and  stimula¬ 
tion  of  Nod2  results  in  the  activation  of  NF-kB  and  MAP  kinase-signaling  cascades  and  production  of  inflam¬ 
matory  molecules  and  anti-microbial  peptides11-13.  Nod2  promotes  a  healthy  microbiota  in  the  intestine,  and 
Nod2- deficient  mice  have  an  altered  gut  microbiota,  which  promotes  colitis14, 15.  Nod2  was  the  first  susceptibility 
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gene  identified  for  Crohn’s  disease16-18,  however,  the  functional  role  of  Nod2  in  the  pathogenesis  of  intestinal 
inflammation  remains  unknown  and  the  role  of  Nod2  in  other  inflammatory  and  non-inflammatory  disorders  is 
mostly  unexplored19-22. 

We  demonstrate  in  this  study  that  deletion  of  the  pattern  recognition  receptor  Nod2  abolishes  the  resistance  of 
BALB/c  to  HFD-induced  obesity  Nod2~'~  BALB/c  mice  on  HFD  gain  substantially  more  weight  than  wild  type 
(WT)  BALB/c  mice  on  HFD  and  become  obese.  Obese  Nod2~'~  BALB/c  mice  develop  metabolic  dysfunction 
including  hyperlipidemia,  hyperglycemia,  and  steatosis,  and  have  increased  adipose  tissue  mass  and  large  lipid 
droplets  in  hepatocytes.  Nod2~'~  mice  on  HFD  have  dysbiosis  and  the  altered  gut  microbiota  contributes  to  the 
development  of  diet- dependent  obesity  and  metabolic  dysfunction. 

Results 

Nod2  protects  from  diet-induced  obesity  in  BALB/c  mice.  WT  and  Nod2~'~  female  mice  on  BALB/c 
background,  matched  for  weight  (average  18g)  and  age  (6-7  weeks  old)  at  the  beginning  of  the  experiment  were 
maintained  on  HFD  (fat  60kcal%)  or  low  fat  diet  (LFD,  fat  10kcal%)  for  30  weeks.  The  total  calories/g  of  food 
was  the  same  for  both  diets.  We  observed  a  dramatic  increase  in  weight  gained  by  Nod2~'~  mice  compared  with 
WT  mice  on  HFD  (Fig.  la).  Nod2~'~  mice  gained  significantly  more  weight  than  WT  mice  starting  from  week  2 
on  HFD,  and  at  week  30  Nod2~'~  mice  gained  almost  twice  as  much  weight  as  WT  mice.  As  expected,  WT  and 
Nod2~'~  mice  on  HFD  gained  significantly  more  weight  than  WT  and  Nod2~'~  mice  on  LFD.  There  was  no  dif¬ 
ference  in  weight  gained  between  WT  and  Nod2~'~  mice  on  LFD.  Our  data  with  WT  mice  concur  with  previous 
data  showing  that  WT  BALB/c  mice  do  not  become  obese  on  HFD4. 

Nod2~'~  mice  on  HFD  were  visibly  larger  with  substantially  more  fat  deposits  than  WT  HFD  mice  (Fig.  lb) 
and  had  significantly  increased  subcutaneous  and  visceral  adipose  tissue  and  an  enlarged  and  fatty  liver  com¬ 
pared  with  WT  HFD  mice  (Fig.  lc  and  d).  There  was  no  difference  in  the  amount  of  food  consumed  or  water 
intake  between  WT  and  Nod2~'~  HFD  mice;  average  food  consumed/mouse/day  was  2.0  g  for  both  genotypes 
and  (Supplementary  Fig.  SI);  the  average  water  intake/mouse/day  was  4.0  ml  for  both  genotypes.  There  was  also 
no  significant  difference  in  the  average  running  distance/24  h  between  WT  (9.48  km)  and  Nod2~/~  (9.59  km)  mice 
on  regular  chow  (Supplementary  Fig.  SI). 

Thus,  our  data  demonstrate  that  in  WT  mice  Nod2  protects  from  the  development  of  HFD-induced  obesity 
and  that  Nod2  deletion  reverses  the  resistance  to  obesity  observed  in  WT  BALB/c  mice.  The  increased  weight 
gain  in  Nod2~'~  HFD  mice  compared  with  WT  HFD  mice  is  not  due  to  increased  food  consumption  or  due  to  an 
inherent  lower  activity 

Nod2  protects  from  diet-induced  hyperlipidemia,  hyperglycemia,  and  glucose  intolerance 
in  BALB/c  mice.  We  next  examined  the  effect  of  Nod2  and  HFD  on  lipid  and  carbohydrate  metabolism. 
Nod2~'~  HFD  mice  exhibited  significantly  higher  levels  of  serum  cholesterol  and  triglycerides  compared  with 
WT  HFD  mice  and  compared  with  Nod2~'~  and  WT  mice  on  LFD  (Fig.  le  and  f).  Nod2~'~  HFD  mice  had  sig¬ 
nificantly  higher  fasting  glucose  and  impaired  ability  to  restore  blood  glucose  in  response  to  a  bolus  of  glucose 
compared  with  WT  HFD  mice  (Fig.  lg  and  h),  which  indicates  development  of  glucose  intolerance.  Nod2~'~  HFD 
mice  had  significantly  higher  levels  of  serum  insulin  than  WT  HFD,  Nod2~'~  LFD,  and  WT  LFD  mice  (Fig.  li), 
and  exhibited  a  mild  but  significantly  reduced  response  to  exogenous  insulin  compared  with  WT  HFD  mice 
120  min  after  an  insulin  challenge  (Fig.  lj),  which  indicates  the  development  of  insulin  resistance  and  suggests 
that  partial  compensation  maybe  achieved  by  increased  insulin  production.  Nod2~'~  HFD  mice  had  significantly 
higher  serum  levels  of  the  anorexic  hormone  leptin  than  WT  HFD  mice,  Nod2~'~  LFD  mice,  and  WT  LFD  mice 
(Fig.  Ik).  WT  HFD  mice  had  significantly  higher  triglycerides  than  Nod2~'~  LFD  and  WT  LFD  mice,  which  indi¬ 
cates  the  development  of  mild  metabolic  dysfunction  in  WT  HFD  mice.  Nod2~'~  LFD  mice  also  had  a  smaller  but 
significant  increase  in  serum  lipids,  glucose,  insulin,  and  leptin  compared  with  WT  LFD  mice.  Thus,  our  results 
demonstrate  that  Nod2~'~  HFD  mice  develop  hyperlipidemia,  hyperglycemia,  and  glucose  intolerance,  which  are 
hallmark  manifestations  of  metabolic  syndrome  in  humans. 

To  determine  any  differences  in  the  total  energy  extracted  from  food  between  Nod2~'~  HFD  mice  and  WT 
HFD  mice,  we  measured  the  total  calories  present  in  the  stool  samples  of  these  mice.  Our  results  show  that  total 
caloric  value  of  Nod2~'~  HFD  feces  was  significantly  lower  than  WT  HFD  feces  (Fig.  11).  These  results  indicate 
that  Nod2~'~  mice  on  HFD  extract  more  energy  from  the  same  amount  of  food  than  WT  mice  on  HFD,  which 
may  be  due  to  differences  in  the  metabolism  of  host  or  gut  microbiota. 

Nod2  protects  from  diet-induced  histopathology  in  BALB/c  mice.  We  characterized  the  histopa- 
thology  in  the  liver  and  adipose  tissue  of  Nod2^/~  obese  mice.  Nod2~/~  HFD  mice  had  significantly  larger  hepat¬ 
ocytes  than  WT  mice  on  HFD  and  the  majority  of  these  cells  had  a  foamy  appearance,  which  indicates  lipid 
accumulation  (Fig.  2a).  In  BODIPY  stained  sections  there  was  a  dramatic  and  significant  increase  in  the  size  of 
lipid  droplets  (LDs)  in  Nod2~'~  mice  compared  with  WT  mice  (Fig.  2b).  The  area  for  LDs  ranged  from  2  \im2  to 
300  [im2  with  the  highest  percentage  of  LDs  in  the  36  to  75  [im2  range  for  Nod2~'~  mice  and  5  to  9  [im2  range  for 
WT  mice  (Fig.  2c).  Nod2~'~  mice  on  HFD  had  significantly  higher  levels  of  liver  cholesterol  and  triglycerides  than 
WT  HFD,  Nod2~'~  LFD,  and  WT  LFD  mice  (Fig.  2d  and  e).  Thus,  our  histological  and  biochemical  data  further 
support  our  results  showing  development  of  obesity  in  Nod2~'~  HFD  mice  and  also  support  our  data  for  WT 
mice,  which  do  not  become  obese  and  do  not  have  a  fatty  liver  (Fig.  1). 

We  next  identified  differences  in  the  proteins  associated  with  LDs  between  Nod2~'~  and  WT  mice  on  HFD. 
We  isolated  proteins  from  LD- enriched  fractions  and  identified  them  by  mass  spectrometry  All  proteins  that 
were  significantly  different  (P<  0.05  and  fold  change  >1.5  or  <0.6)  between  Nod2~'~  and  WT  mice  are  directly 
or  indirectly  involved  in  lipid  metabolism  or  in  vesicle  trafficking  (Fig.  2f).  Nod2~'~  HFD  mice  compared  with 
WT  HFD  mice  had  significantly  higher  levels  of  Abhd5,  Agpat2,  Gpam,  and  Gpd2  but  decreased  levels  of  Gpat4, 
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Figure  1.  Nod2~'~  mice  on  HFD  become  obese  and  develop  hyperlipidemia,  hyperglycemia,  and  glucose 
intolerance.  WT  and  Nod2~'~  mice  were  monitored  for  weight  gain,  increased  adiposity,  and  assayed  for 
metabolites,  (a)  Percent  gain  in  body  weight  over  week  0.  (b-d)  At  week  30,  (b)  representative  mice  are  shown 
for  increase  in  overall  size  and  visceral  adipose  tissue,  (c)  average  increase  in  subcutaneous  thickness  and  weight 
of  visceral  adipose  tissue  and  liver  per  mouse,  and  (d)  total  visceral  adipose  tissue  and  liver  from  representative 
mice,  (e-j)  At  week  30,  mice  were  assayed  for  serum  (e)  cholesterol,  (f)  triglycerides,  (g)  glucose,  (i)  insulin,  and 
(k)  leptin,  and  for  the  development  of  (h)  glucose  tolerance  and  (j)  insulin  resistance.  (1)  Total  caloric  value  in 
fecal  samples.  The  results  are  (a)  means  ±  SEM  of  18-25  mice/group;  (c,h,j)  means  ±  SEM  of  10-12  mice/group; 
or  (e-g,i,k,l)  for  individual  mice/group  ±  SEM.  (h)  Data  for  the  glucose  tolerance  test  is  shown  as  ratio  of  HFD 
to  LFD  for  each  strain.  *P<  0.05,  **P<  0.001,  Nod2~'~  HFD  versus  WT  HFD;  *P<  0.05,  ##P<  0.001,  Nod2~/~ 
HFD  versus  Nod2~l~  LFD;  AP  <  0.05,  AAP<  0.001,  WT  HFD  versus  WT  LFD;  +P<  0.05,  ++P<  0.001,  Nod2~'- 
LFD  versus  WT  LFD.  The  results  for  food  consumed  and  voluntary  wheel  activity  are  shown  in  Fig.  SI. 
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Figure  2.  Nod2~'~  mice  on  HFD  develop  steatosis  and  have  large  stressed  adipocytes  with  increased  infiltration 
of  macrophages.  (a,b)  Liver  sections  from  Nod2~'~  and  WT  mice  stained  with  (a)  H&E  or  (b)  BODIPY  and 
DAPI.  (c)  Percent  frequency  distribution  of  the  area  (pm2)  of  -5,000  LDs/group.  (d,e)  Quantification  of  liver 
cholesterol  and  triglycerides,  (f)  Heatmap  representation  of  the  fold  increase  (red)  or  decrease  (green)  in 
the  abundance  of  proteins  present  in  LD -enriched  fraction  from  liver.  The  fold  ratio  for  individual  Nod2~'~ 

HFD  mice  to  the  average  of  WT  HFD  mice  were  computed  (lanes  1  to  5)  with  the  average  of  the  fold  for  all 
Nod2~'~  HFD  mice,  (g)  Adipose  tissue  H&E  sections  and  a  dying  cell  (— >)  is  indicated,  (h)  Percent  frequency 
distribution  the  area  of -8,000  adipocytes/group,  (i)  Average  numbers  of  dying  adipocytes  (crowns)  per 
3,000  cells  in  H&E  stained  sections,  (j)  Quantification  of  CD45+F4/80+  macrophage  in  adipose  tissue,  (k-m) 
Transcript  levels  for  Adgrel ,  Cd68 ,  and  Tnfa  in  adipose  tissue.  Images  are  representative  and  results  are  (c,h,i) 
means  ±  SEM  of  6  mice/group  and  (d,e,j-m)  for  individual  mice  with  N=  5-10  mice/group.  *P  <  0.05, 

**P<  0.001,  Nod2~'~  HFD  versus  WT  HFD;  *P<  0.05,  ##P<  0.001,  Nod2~'~  HFD  versus  Nod2~l~  LFD; 

AP  <  0.05,  WT  HFD  versus  WT  LFD;  and  +P  <  0.05,  Nod2~l~  LFD  versus  WT  LFD. 
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proteins  associated  with  biosynthesis  of  glycerolipids.  Tecr,  a  protein  involved  in  the  synthesis  of  very  long  chain 
fatty  acids  was  increased  in  Nod2~'~  HFD  mice.  Many  proteins  that  participate  in  the  catabolism  of  lipids  were 
also  more  abundant  in  Nod2~'~  HFD  mice  than  in  WT  HFD  mice,  including  Acaalb,  Acatl,  Aldh7al,  Bdhl, 
Cyp2b9,  Cyp4al4,  Faah,  Lipe,  Lyplal,  Mgll,  and  Pnpla3.  Proteins  associated  with  lipid  storage  were  selectively 
up  regulated  or  down  regulated,  Plin3  and  Plin4  were  present  at  higher  and  Faf2  and  Plinl  at  lower  levels  in 
Nod2~'~  HFD  mice.  Several  proteins  involved  in  vesicle  trafficking  were  decreased  in  the  LD-enriched  fraction 
of  Nod2~/~  HFD  mice  compared  with  WT  HFD  mice.  However,  Copbl,  a  protein  associated  with  non-clathrin 
coated  vesicles  with  a  role  in  limiting  lipid  storage,  was  significantly  higher  in  Nod2~'~  HFD  mice  than  WT  HFD 
mice.  These  data  indicate  that  the  larger  LDs  in  Nod2~'~  HFD  mice  compared  with  WT  HFD  mice  are  accompa¬ 
nied  by  changes  in  lipid  metabolizing  enzymes  and  in  vesicle  trafficking  proteins. 

Nod2~'~  mice  on  HFD  had  enlarged  adipocytes  with  significantly  larger  area  than  WT  HFD  mice  (Fig.  2g  and 
h).  The  area  for  adipocytes  ranged  from  31  |im2  to  16,000  |im2  with  the  highest  percentage  of  cells  in  the  2000  to 
4000  [im2  range  for  Nod2~'~  mice  and  250  to  500  [im2  for  WT  mice  (Fig.  2h).  Nod2~'~  mice  had  significantly  more 
“crowns”,  which  are  dead  adipocytes  surrounded  by  immune  cells  (Fig.  2g  and  i,  indicated  by  an  arrow)  and  are  a 
characteristic  feature  of  stressed  adipose  tissue  associated  with  obesity  Nod2~'~  mice  on  HFD  had  significantly 
more  F4/80+  macrophages  in  adipose  tissue  than  WT  mice  on  HFD,  Nod2~'~  LFD,  and  WT  LFD  mice  (Fig.  2j) 
with  a  corresponding  increase  in  the  transcripts  for  Adgrel  (gene  for  F4/80)  and  Cd68 ,  another  marker  for  mono¬ 
cytes/macrophages  (Fig.  2k  and  1).  Expression  of  Tnfa ,  an  inflammatory  cytokine  associated  with  obesity  was 
enhanced  (Fig.  2m).  Thus,  our  data  demonstrate  that  Nod2  protects  from  lipid  accumulation  in  the  adipose  tissue 
and  liver,  and  that  development  of  diet-dependent  obesity  in  Nod2~'~  HFD  mice  is  associated  with  steatosis, 
formation  of  large  LDs  in  hepatocytes,  stressed  adipocytes,  and  infiltration  of  macrophages  in  the  adipose  tissue. 

Diet-dependent  obesity  in  Nod2~'~  mice  is  associated  with  changes  in  the  expression  of  adi¬ 
pose  tissue  genes  involved  in  immune  responses  and  intermediary  metabolism.  We  identified 
genes  that  were  differentially  expressed  in  the  visceral  adipose  tissue  of  WT  and  Nod2~'~  BALB/c  mice  main¬ 
tained  on  HFD  for  30  weeks  using  RNA  transcriptomics.  Nod2~/~  HFD  mice  had  357  genes  with  significantly 
increased  expression  (and  fold  change  >2)  and  169  genes  with  significantly  decreased  expression  (and  fold 
change  <0.5)  compared  with  WT  mice  on  HFD  (Supplementary  Fig.  S2).  79%  of  the  differentially  expressed 
genes  are  associated  with  a  function  or  pathway  (Ensembl  or  DAVID)  and  are  shown  in  Fig.  3  and  Supplementary 
Table  SI.  More  than  100  genes  that  directly  or  indirectly  participate  in  immune  responses  were  differentially  reg¬ 
ulated  and  the  majority  of  these  genes  (>80)  were  significantly  up  regulated  in  the  adipose  tissue  of  Nod2~'~  HFD 
mice  compared  with  WT  HFD  mice.  These  genes  code  for  transcription  factors,  signaling  molecules,  effector 
molecules,  enzymes,  and  receptors,  all  involved  in  regulating  immune  cells  and  responses.  For  example,  Nod2~'~ 
mice  on  HFD  had  increased  expression  of  genes  for  (i)  chemokines  Cc/2,  Cc/5,  Cc/7,  Ccl8,  Ccll2 ,  Cc/22,  and  Cxcll , 
which  recruit  monocytes,  eosinophils,  and  lymphocytes;  (ii)  enzymes  involved  in  arachidonic  acid  and  eicosa- 
noid  metabolism,  Cyp2c39 ,  Cyp2dl0,  Cyp4fl4 ,  Fam213b ,  Pla2gl2b,  and  Ptges3l ;  (iii)  receptors  involved  in  T  cell 
activation,  including  Cd36 ,  Cd3e ,  Cd3y,  Cd4 ,  Cd5,  Cd6,  and  Cd52 ;  and  (iv)  adipokines  Sfrp5,  Stra6l ,  and  Wisp2 , 
and  the  protease  Cpnl.  Expression  of  genes  for  adipsin  ( Cfd )  and  regulators  of  adipokines,  Hcar2  and  Ptprn2 ,  was 
decreased. 

The  second  largest  group  of  genes  that  was  differentially  regulated  between  WT  and  Nod2~'~  HFD  mice 
directly  or  indirectly  participates  in  intermediary  metabolism,  primarily  lipid  metabolism  (Fig.  3,  Supplementary 
Table  SI).  These  genes  include:  (i)  enzymes  involved  in  breakdown  of  fatty  acids  and  glycerolipids,  for  example, 
Acot3 ,  Acot4 ,  Ces4a ,  Agpat9 ,  and  Pnplal  were  increased;  (ii)  members  of  the  acyl-CoA  synthetase  medium  chain 
family  were  up  regulated  ( Acsml )  or  down  regulated  ( Acsm3  and  Acsm5);  (iii)  enzymes  in  the  synthesis  of  fatty 
acids  and  triglycerides,  Mogat2  and  Scdl ,  were  increased;  (iv)  transporters  for  fatty  acids,  cholesterol,  glycerolip¬ 
ids,  and  glycerol,  Abcc2 ,  Aqpll ,  Atp8b5 ,  Fabp5 ,  and  Mfge8 ,  were  up  regulated,  and  Apoldl ,  Glccil ,  and  Slc27al , 
were  down  regulated. 

Nod2~'~  mice  on  HFD  had  several  differentially  expressed  genes  that  directly  or  indirectly  participate  in 
amino  acid  and  carbohydrate  metabolism  compared  with  WT  HFD  mice  (Fig.  3,  Supplementary  Table  SI).  These 
genes  include:  (i)  enzymes  in  amino  acid  metabolism,  Accsl ,  Agxt ,  and  Aspg  were  decreased  or  Got  111,  Hgd,  Sds, 
and  Sdsl  were  increased  and  (ii)  amino  acid  transport,  Slc22a2 ,  Slc7al4 ,  and  Urocl  were  increased.  Genes  for:  (i) 
enzymes  in  carbohydrate  metabolism,  Gys2 ,  Lctl ,  Mgam ,  and  Ppplr3e  and  (ii)  glucose  transport,  AdamtsS,  Rtn2, 
Slc2a2 ,  Smpd5,  and  Syp  were  up  regulated,  but  (iii)  insulin  signaling,  Dusp9,  Fbxo40 ,  Grbl4 ,  and  Ptpre  were  down 
regulated. 

Nod2~'~  mice  on  HFD  had  an  overall  decrease  in  the  expression  of  genes  involved  in  (i)  neural  signaling, 
including  Brinp3 ,  Cadps2 ,  Dlg2 ,  Itpka ,  Mylka ,  Prickle2 ,  Shisa6,  Sores  1,  Sytl5 ;  (2)  nerve  adhesion,  Celsr2,  Cntnl, 
Fat3 ,  Lsamp ,  Ncaml ,  Ncam2 ,  Nfasc,  and  Nlgnl;  and  (iii)  transcriptional  regulators  for  neural  development,  Ascii, 
Hesxl,  and  Id4  (Fig.  3,  Supplementary  Table  SI). 

Genes  involved  in  other  cellular  functions,  including  cell  cycle,  development,  cell  adhesion,  and  cytoskele- 
ton  were  also  differentially  expressed  (Fig.  3,  Supplementary  Table  SI).  Nod2~'~  HFD  mice  also  had  an  over¬ 
all  decrease  in  the  expression  of  genes  involved  in  (i)  cell  adhesion  and  signaling,  including  Adgral,  Adgrg2, 
Cdhl2,  Cdh9,  Cldn22,  Coll2al,  Col6a5,  Col6a6,  Pcdh20,  Pcdhbl2,  Pcdhbl6,  Pcdhbl8,  Pcdhb9,  Rassf6,  Snedl ,  and 
Thsd4;  and  (ii)  inhibition  of  apoptosis,  including  Angptl,  Dusp26,  Egln3,  Hspala ,  and  Pde3a.  There  was  an  over¬ 
all  increase  in  expression  of  genes  involved  in:  (i)  DNA  structure,  repair  and  chromatin  silencing,  for  exam¬ 
ple,  Histlhld,  Histlhle,  Histlh2al,  Histlh2be ,  Hist2h3c2,  and  Mutyh;  (ii)  transcription  and  post-transcriptional 
processing,  Abra,  Inpp4b,  Lbx2,  Pcbdl,  Pnldcl ,  Rbm28,  and  Rpphl;  and  (iii)  proteases  and  protease  inhibitors, 
Adam23,  Cpa2,  Itihl,  Itih3,  Lonrf3,  Psma8,  Serpinall ,  Serpinald ,  Serpinale,  Serpina3e-ps,  Serpina3m,  and 
Wfdcl8.  However,  there  was  a  noticeable  decrease  in  genes  involved  in  translation  and  post-translational  pro¬ 
cessing,  B3galtl,  B3galt2,  Faml29a,  GalntS,  St8sia2,  Tgm3,  TgmS,  Ttll2,  and  Ttll7. 
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Figure  3.  Differential  expression  of  genes  for  immunity,  metabolism,  and  other  cellular  functions  in  the 
adipose  tissue  of  Nod2~'~  mice  on  HFD.  Heatmap  representation  of  the  fold  increase  (red)  or  decrease  (green) 
for  gene  expression  in  the  adipose  tissue  of  Nod2~/~  HFD  mice  compared  with  WT  HFD  mice.  The  genes  are 
grouped  based  on  function  or  pathway  The  fold  ratio  for  individual  Nod2~'~  HFD  mice  to  the  average  of  WT 
HFD  mice  was  computed  and  is  shown  (lanes  1  to  6)  with  the  average  of  the  fold  for  all  Nod2~'~  HFD  mice. 
Genes  that  had  a  fold  change  of  >2  or  <0.5  (<— 2)  and  P  <  0.05  with  5%  FDR  are  included  in  the  heatmap. 
N=  6  mice/group.  The  numerical  data  for  fold  increase  in  individual  mice,  the  average,  P  value,  and  FDR  are 
shown  in  Supplementary  Table  SI. 
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Thus,  our  results  indicate  that  genes  involved  in  a  wide  array  of  pathways  and  cellular  functions  are  differen¬ 
tially  regulated  in  the  adipose  tissue  of  Nod2~'~  HFD  mice  compared  with  WT  HFD  mice.  These  changes  in  gene 
expression  may  contribute  to  the  increased  adiposity,  steatosis,  obesity,  and  metabolic  dysfunction  observed  in 
Nod2~'~  mice  on  HFD.  Future  studies  should  confirm  these  changes  at  the  protein  level  and  focus  on  their  role  in 
the  development  of  obesity  in  Nod2~'~  HFD  mice. 

Diet-dependent  obesity  in  Nod2~'~  mice  is  associated  with  changes  in  the  expression  of  liver 
genes  involved  in  intermediary  metabolism.  We  next  identified  genes  that  were  differentially  expressed 
in  the  liver  of  WT  and  Nod2~'~  BALB/c  mice  maintained  on  HFD  for  30  weeks  using  RNA  transcriptomics. 
Nod2~'~  HFD  mice  had  188  genes  with  significantly  increased  expression  (and  fold  change  >2)  and  130  genes 
with  significantly  decreased  expression  (and  fold  change  <0.5)  compared  with  WT  mice  on  HFD  (Supplementary 
Fig.  S2).  87%  of  the  differentially  expressed  genes  are  associated  with  a  function  or  pathway  (Ensembl  or  DAVID) 
and  are  shown  in  Fig.  4  and  Supplementary  Table  S2.  Of  these  genes,  there  were  more  than  60  genes  that  directly 
or  indirectly  participate  in  intermediary  metabolism,  primarily  lipid  metabolism  and  include  genes  for  enzymes, 
transport,  and  regulation.  For  example,  Nod2~'~  HFD  mice  had  higher  expression  of  enzymes  and  regulators 
involved  in:  (i)  synthesis  and  storage  of  fatty  acids,  triglycerides,  steroids,  and  glycolipids,  Acacb ,  Acss2y  Elovl3 , 
Elovl6 ,  Fitml ,  G0s2 ,  Hsd3b5 ,  Mogatl ,  Scdl ,  and  St6galnac4 ;  (ii)  catabolism  of  lipids,  Pla2g4f  Pnpla3 ,  Pnpla5 , 
Pexlla ,  and  Phyhip ;  and  (iii)  carbohydrate  metabolism,  Gck ,  Pdk4 ,  and  Ppplr3e.  Several  genes  involved  in  trans¬ 
port  of  metabolites  and  nutrients  were  down  regulated,  including  Aqp4,  Fabp5 ,  Ndgrl ,  and  SlcolAl,  whereas 
other  transporters  were  up  regulated,  including  Cd36 ,  Mfsd2a>  Osbpl5 ,  Pltp ,  Slcl0a2 ,  Slcl6al3 ,  and  Slc2a4  (Glut4). 
Nod2~'~  HFD  mice  had  decreased  expression  of  the  leptin  receptor  ( Lepr ). 

Nod2~'~  mice  on  HFD  had  about  50  genes  involved  in  immune  responses  that  were  differentially  expressed 
in  the  liver  compared  with  WT  mice  on  HFD  (Fig.  4,  Supplementary  Table  S2).  Transcription  factors  involved  in 
the  activation  of  immune  cells,  Arl2bp ,  Cbfa2t3 ,  and  Cebpe ,  were  increased,  whereas  transcriptional  repressors 
Bcl6  and  Camlkd  were  decreased.  There  was  a  noticeable  increase  in  the  expression  of  many  interferon-induced 
genes,  including  the  chemokines  Cxc/9,  CxcllO ,  and  Cxclll.  Genes  for  the  T  cell  receptor  Cd8a  and  retinoic  acid 
metabolism,  Cyp26bl,Raetl6 ,  and  Raetlz ,  that  regulate  T  and  NK  cell  immune  responses  were  induced,  but  genes 
in  arachidonic  acid  metabolism,  Aloxl5 ,  Cyp2c54 ,  Cyp2gl ,  and  Cyp2sl ,  were  inhibited. 

Nod2~/~  mice  on  HFD  had  increased  expression  of  many  liver  genes  involved  in  DNA  structure,  replication, 
and  repair  (Fig.  4,  Supplementary  Table  S2).  Genes  for  serine  protease  inhibitors  ( ItihS ,  Serpina4-psl ,  Serpina7 , 
and  Serpina9 ),  chaperones  ( Hspla  and  Hsplb ),  and  ubiquitination  ( Kelch33 ,  Klhl33 ,  Triml3 ,  and  Ube2qll)  were 
decreased.  Genes  for  cell  cycle  and  cell  proliferation  were  overwhelmingly  induced,  however  many  genes  involved 
in  development  were  down  regulated,  including  transcription  factors  (Jrxl,  K7/4,  Meoxl ,  Sox9 ,  Sp5 ,  and  Tbx3 ), 
signaling  molecules  ( Arhgef40 ,  Dact2 ,  Mpzll ,  Mrgpre ,  and  Pdgfc)  and  components  of  the  Wnt  pathway  ( Axin2 , 
Bcl7c ,  Wifi,  and  Lgr5).  Several  genes  involved  in  neural  development  were  down  regulated,  including  proteins 
involved  in  ion  and  vesicular  transport  ( Adamll ,  Cacnald ,  Syngrl ,  and  Syf9)  and  neural  signaling  ( Aphlc , 
Cadm4 ,  and  Pde6c).  However,  genes  for  acetylcholine  receptor  ( Chrna4 )  and  signaling  at  the  neuronal  junction 
were  up  regulated  ( Clstn3 ,  Cntnapl ,  and  Pclo).  The  tyrosine  kinase  receptor,  Ntrk2 ,  which  has  a  role  in  suppress¬ 
ing  appetite,  was  also  up  regulated  (Fig.  4,  Supplementary  Table  S2). 

Nod2~'~  mice  on  HFD  had  increased  expression  of  many  genes  involved  in  extracellular  matrix  organization 
and  cell  adhesion,  including  Col27al ,  Esml ,  Haplnl ,  Hmcn2 ,  Lamb3 ,  Myo7b ,  and  Pcdhl7.  Genes  involved  in  ves¬ 
icle  trafficking  ( Rabl  lfip4 ,  Rab2b ,  Rasl2-9 ,  and  Spffr)  and  proton,  sodium  and  sulfur  transport  ( Atp6v0d2 ,  Slc4a9 , 
Slc8a3 ,  Sult3al ,  Sult5al ,  and  Wnk4)  were  down  regulated,  whereas  genes  for  zinc  and  iron  uptake  ( Slc39a5  and 
Tfrc)  were  increased  (Fig.  4,  Supplementary  Table  S2). 

These  results  indicate  that  in  the  liver,  genes  involved  in  a  wide  array  of  pathways  and  cellular  functions  are 
differentially  regulated  by  HFD  in  Nod2~'~  mice  compared  with  WT  mice.  These  changes  in  gene  expression  may 
contribute  to  the  increased  adiposity,  steatosis,  obesity,  and  metabolic  dysfunction  observed  in  Nod2~'~  mice  on 
HFD.  Future  studies  should  confirm  these  changes  at  the  protein  level  and  focus  on  their  role  in  the  development 
of  obesity  in  Nod2~'~  HFD  mice. 

Nod2~'~  genotype  and  HFD  and  LFD  affect  diversity  of  stool  microbiome.  Because  intestinal 
microbiome  is  important  for  energy  extraction  from  food,  we  determined  the  diversity  of  intestinal  bacteria  in 
WT  and  Nod2~'~  mice  on  HFD  and  LFD,  to  identify  in  detail  the  differences  in  their  intestinal  microbiota.  We 
isolated  DNA  from  stool  microbiota  and  performed  genetic  phylotyping  (community  profiling)  using  pyrose- 
quencing  of  the  variable  regions  of  bacterial  16S  ribosomal  RNA  (rRNA)  genes.  To  determine  the  effect  of  geno¬ 
type,  we  compared  microbiomes  in  WT  and  Nod2~'~  mice  on  the  same  diets,  and  to  determine  the  effect  of  diets, 
we  compared  microbiomes  in  mice  of  the  same  genotype  fed  HFD  or  LFD. 

We  first  evaluated  a-diversity  of  the  microbiomes  of  WT  and  Nod2~'~  mice  on  HFD  and  LFD.  We  identified 
a  total  of  223  species  and  4425  operational  taxonomic  units  (OTUs)  in  the  stools  of  all  mice:  153,  172,  180,  and 
164  species  and  2344, 2589, 2392,  and  2381  OTUs  in  WT  HFD,  Nod2~'~  HFD,  WT  LFD,  and  Nod2~'~  LFD  mice, 
respectively  The  total  number  of  OTUs  identified  in  Nod2~'~  HFD  mice  was  significantly  higher  than  in  WT  HFD 
mice,  whereas  the  total  numbers  of  OTUs  in  Nod2~'~  LFD  and  WT  LFD  mice,  as  well  as  total  numbers  of  species 
in  all  groups  of  mice  were  similar  with  no  significant  differences  between  the  groups  (Supplementary  Fig.  S3).  The 
numbers  of  species  and  OTUs  per  mouse  (microbiome  richness)  were  similar  in  all  groups  of  mice,  except  for  WT 
LFD  mice,  which  had  significantly  more  species/mouse  than  WT  HFD  mice.  Both  Shannon  diversity  indices  (H) 
and  Shannon  equitability  indices  (EH,  reflecting  microbiome  evenness)  were  significantly  higher  for  the  species 
in  Nod2~'~  HFD  than  in  WT  HFD  mice,  and  for  both  the  species  and  OTUs  in  WT  LFD  than  in  Nod2~ LFD 
mice,  and  also  for  both  the  species  and  OTUs  in  WT  LFD  than  in  WT  HFD  mice  (Supplementary  Fig.  S3).  These 
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Figure  4.  Differential  expression  of  genes  for  immunity,  metabolism,  and  other  cellular  functions  in  the  liver 
of  Nod2~'~  mice  on  HFD.  Heatmap  representation  of  the  fold  increase  (red)  or  decrease  (green)  for  gene 
expression  in  the  adipose  tissue  of  Nod2~'~  HFD  mice  compared  with  WT  HFD  mice.  The  genes  are  grouped 
based  on  function  or  pathway  The  fold  ratio  for  individual  Nod2~'~  HFD  mice  to  the  average  of  WT  HFD  mice 
was  computed  and  is  shown  (lanes  1  to  6)  with  the  average  of  the  fold  for  all  Nod2~'~  HFD  mice.  Genes  that 
had  a  fold  change  of>2or<0.5(<  — 2)  and  P  <  0.05  with  5%  FDR  are  included  in  the  heatmap.  N=  6  mice/ 
group.  The  numerical  data  for  fold  increase  in  individual  mice,  the  average,  P  value,  and  FDR  are  shown  in 
Supplementary  Table  S2. 
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results  show  that  both  Nod2~'~  genotype  on  HFD  and  WT  genotype  on  LFD  may  increase  a-diversity,  although 
these  changes  are  very  modest. 

We  then  evaluated  (3-diversity  of  the  microbiomes.  Microbiota  from  stools  of  all  four  groups  of  mice  (WT 
HFD,  Nod2~'~  HFD,  WT  LFD  and  Nod2~'~  LFD)  separated  from  each  other  in  the  Principal  Coordinate  Analysis 
(PCoA),  with  microbiota  from  each  group  of  mice  showing  distinct  and  significant  separation  from  the  microbi¬ 
ota  of  all  other  groups  of  mice  (Fig.  5a).  These  results  suggested  differences  in  (3-diversity  in  the  intestinal  micro¬ 
biomes  between  all  four  groups  of  mice  with  both  the  genotype  (WT  and  Nod2~'~)  and  the  diet  (HFD  and  LFD) 
affecting  (3-diversity  of  stool  microbiomes.  These  results  prompted  us  to  further  identify  the  differences  between 
these  microbiomes. 

We  detected  many  statistically  significant  effects  of  both  Nod2  genotype  and  the  diet  on  bacterial  (3-diversity, 
with  many  significant  changes  in  the  abundance  of  various  bacterial  groups  in  the  stools  at  all  taxonomic  levels. 
For  example,  at  the  class  level  (Fig.  5b  and  Supplementary  Fig.  S4),  WT  HFD  mice  had  significantly  increased 
abundance  of  Epsilonproteobacteria  compared  with  Nod2~'~  HFD  mice,  and  Nod2~'~  HFD  mice  had  significantly 
increased  abundance  of  Deferribacteres  compared  with  WT  HFD  mice.  In  mice  on  LFD,  WT  mice  had  signifi¬ 
cantly  higher  abundance  of  Candidatus  saccharibacteria  than  Nod2~'~  mice,  and  Nod2~'~  mice  had  significantly 
higher  abundance  of  Deferribacteres  than  WT  mice.  The  diet  also  had  a  pronounced  effect  on  stool  microbiota. 
HFD  significantly  increased  the  abundance  of  Clostridia  and  Deltaproteobacteria  in  both  WT  and  Nod2~'~  mice 
and  also  Deferribacteres  in  Nod2~'~  mice,  whereas  LFD  significantly  increased  the  abundance  of  Bacteroida , 
Bacilli ,  and  Erysipelotrichia  in  both  WT  and  Nod2~'~  mice,  and  also  Candidatus  saccharibacteria  in  WT  mice. 

At  other  taxonomic  levels,  significant  differences  in  abundance  in  the  following  numbers  of  orders,  families, 
genera,  species,  and  OTUs  (out  of  total  of  27,  47,  109,  223,  and  4425)  were  detected:  between  WT  HFD  and 
Nod2~'~  HFD  mice,  2, 2, 9, 26,  and  548;  between  WT  LFD  and  Nod2~'~  LFD  mice,  6, 10, 20, 46,  and  601;  between 
WT  HFD  and  WT  LFD,  8,  14,  24,  48,  and  535;  and  between  Nod2~,~  HFD  and  Nod2~'~  LFD  mice,  6,  9,  16,  35, 
and  288,  respectively  (NCBI  SRA,  accession  No.  SRP076031).  Thus,  both  Nod2  genotype  of  mice  and  their  diet 
together  significantly  influence  the  diversity  of  their  microbiomes. 

We  then  further  analyzed  the  differences  in  microbiomes  between  Nod2~'~  HFD  mice  and  the  remain¬ 
ing  groups  of  mice,  because  these  mice  had  significantly  higher  weight  gain  and  significant  manifestations  of 
metabolic  syndrome  compared  with  the  other  groups  of  mice,  and  moreover,  these  differences  were  abolished 
by  treatment  with  antibiotics  (next  section).  Nod2~'~  HFD  mice,  compared  with  WT  HFD  mice,  had  signif¬ 
icantly  higher  abundance  of  six  Firmicutes ,  three  Bacteroidetes ,  two  Deferribacteres ,  one  Proteobacteria ,  and 
one  Actinobacteria  species  (Table  1).  Importantly,  three  of  these  bacterial  species,  both  Deferribacteres  and  one 
Firmicutes  ( Lachnoclostridium  phytofermentans)  had  significantly  higher  abundance  in  Nod2~'~  HFD  mice  than 
in  the  other  three  groups  of  mice,  thus  showing  the  best  correlation  with  obesity  and  metabolic  syndrome,  since 
Nod2~'~  HFD  mice  was  the  only  group  with  significantly  higher  obesity  and  metabolic  syndrome.  There  were  no 
bacterial  species  with  significantly  decreased  abundance  in  Nod2~'~  HFD  mice  compared  with  all  other  groups 
of  mice. 

Diet-induced  obesity  in  Nod2~'~  BALB/c  mice  is  abolished  by  long-term  antibiotic  treat¬ 
ment.  We  next  determined  the  role  of  the  gut  microbiota  in  diet-dependent  obesity  in  Nod2~'~  mice  by 
depleting  gut  bacteria  with  antibiotics.  WT  and  Nod2~'~  mice,  starting  at  4  weeks  of  age,  were  given  an  antibiotic 
mix  (Abx)  of  ciprofloxacin  and  metronidazole  in  their  drinking  water,  which  is  an  established  method  for  deplet¬ 
ing  gut  bacteria20,23-25.  After  three  weeks  the  animals  were  placed  on  HFD  and  Abx  were  continued,  and  body 
weight  was  measured  each  week  (Fig.  6a).  Control  mice  were  started  on  HFD  at  the  same  age,  but  never  treated 
with  Abx.  Our  results  demonstrate  that  Abx  treatment  completely  abolished  the  diet- dependent  weight  gain  in 
Nod2~'~  mice,  but  did  not  affect  the  weight  gain  in  WT  mice  (Fig.  6b).  Nod2~'~  Abx  HFD  mice  gained  signifi¬ 
cantly  less  weight  than  Nod2~'~  HFD  mice,  and  also  less  than  WT  HFD  mice  with  and  without  Abx.  There  was 
no  difference  in  weight  gained  between  WT  Abx  HFD  and  WT  HFD  mice.  There  was  no  difference  in  food  con¬ 
sumption  or  water  intake  between  the  WT  Abx  HFD  and  Nod2~'~  Abx  HFD,  as  average  food  consumed/mouse/ 
day  was  1.86  g  and  1.89  g,  respectively,  and  the  average  water  intake/ mouse/ day  was  2.9  ml  for  both  genotypes. 
There  is  some  decrease  in  water  intake  between  the  antibiotic  treated  versus  control  mice,  however  this  decrease 
does  not  account  for  the  difference  in  the  weight  gain  observed  between  WT  and  Nod2~'~  mice,  because  both 
genotypes  on  antibiotics  had  similar  decrease  in  water  intake  but  only  Nod2~/~  mice  had  significantly  reduced 
weight  gain. 

We  next  determined  whether  the  antibiotic-induced  loss  in  weight  gain  was  accompanied  by  decreased  lipid 
accumulation  in  the  liver  of  Nod2~'~  mice  using  histological  and  biochemical  assays.  Our  results  demonstrate  that 
Nod2~'~  mice  treated  with  long-term  Abx  and  maintained  on  HFD  had  smaller  LDs  in  hepatocytes  compared 
with  Nod2~'~  mice  not  treated  with  Abx  but  maintained  on  HFD  (Fig.  6c).  The  area  for  LDs  ranged  from  2  pm2 
to  300  pm2  with  the  highest  percentage  of  LDs  in  the  36  to  75  pm2  range  for  the  HFD  mice  and  2  to  19  pm2  range 
for  Abx  +  HFD  mice  (Fig.  6d).  Nod2~'~  mice  on  Abx  +  HFD  had  significantly  lower  levels  of  liver  cholesterol  and 
triglycerides  than  control  Nod2~'~  mice  on  HFD  alone  (Fig.  6e  and  f).  Furthermore,  Nod2~'~  mice  on  Abx  +  HFD 
had  significantly  lower  levels  of  serum  cholesterol  and  triglycerides  (Fig.  6g  and  h),  but  no  significant  difference 
in  glucose,  compared  with  control  Nod2~'~  mice  on  HFD  alone. 

Thus,  our  results  demonstrate  that  the  gut  microbiota  in  Nod2~'~  HFD  mice  is  required  for  the  development 
of  diet-dependent  obesity,  steatosis,  and  hyperlipidemia  in  Nod2~'~  mice,  and  that  long-term  treatment  with 
broad-spectrum  antibiotics  abolishes  weight  gain  and  several  aspects  of  metabolic  dysfunction. 

Intestinal  microbiota  from  Nod2~'~  mice  on  high  fat  diet  increases  sensitivity  to  obesity.  We 

next  considered  whether  the  altered  gut  microbiota  in  Nod2~'~  mice  on  HFD  is  sufficient  for  the  development 
of  diet  dependent  obesity.  We  colonized  germ-free  pregnant  female  WT  Swiss-Webster  mice  by  co-housing  with 
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Figure  5.  (3-diversity  in  stool  microbiota  of  WT  HFD,  Nod2~'~  HFD,  WT  LFD,  and  Nod2^'~  LFD  mice,  (a) 
Principal  Coordinate  Analysis  (PCoA)  by  UniFrac  (unweighted)  of  microbiomes.  Each  sphere  with  confidence 
elipsoid  corresponds  to  a  stool  microbiome  from  one  mouse.  N=  6  mice/group,  (b)  Class  abundance  expressed 
as  %  of  total  stool  microbiota;  *  classes  with  significantly  (at  P<  0.05)  increased  abundance  in  WT  versus 
Nod2~'~  mice  (on  HFD  or  LFD);  #classes  with  significantly  (at  P<  0.05)  increased  abundance  in  HFD  versus 
LFD  groups  (for  WT  or  Nod2~'~  mice);  N=  6  mice/group.  The  results  for  individual  mice  are  shown  in 
Supplementary  Fig.  S4.  The  entire  microbiome  analysis  data  comparing  bacterial  diversity  in  the  stools  of  WT 
and  Nod2~'~  mice  maintained  on  HFD  and  LFD  have  been  deposited  in  NCBI  SRA,  accession  No.  SRP076031. 


either  Nod2~'~  HFD  mice  or  WT  HFD  mice,  and  also  by  weekly  gavaging  with  microbiota  from  either  Nod2~'~ 
HFD  mice  or  WT  HFD  mice.  We  placed  these  mice  on  HFD  after  the  first  gavage  (Fig.  7a).  We  removed  the 
co-housed  mice  one  day  after  delivery,  but  continued  weekly  gavaging  and  HFD,  and  monitored  the  weight  gain 
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Phylum 

Species 

Nod2~'~  HFD  (% 
abundance) 

WT  HFD  (% 
abundance) 

Firmicutes 

Hydrogenoanaerobacterium  sp. 

0.76 

0.52 

Coprococcus  eutactus 

0.34 

0.052 

Eubacterium  rectale 

0.18 

0.049 

Lachnoclostridium  (Clostridium) 
phytofermentans 

0.15 

0.0000 

Coprococcus  catus 

0.076 

0.0000 

Blautia  (Ruminococcus)  obeum 

0.027 

0.0025 

Bacteroidetes 

Rikenella  sp. 

4.49 

1.49 

Bacteroides  sp. 

1.58 

0.71 

Butyricimonas  sp. 

0.011 

0.0000 

Deferribacteres 

Mucispirillum  sp. 

2.27 

0.0000 

Mucispirillum  (Flexistipes)  sp. 

0.30 

0.0000 

Proteobacteria 

Helicobacter  ganmani 

0.058 

0.0000 

Actinobacteria 

Adlercreutzia  equolifaciens 

0.028 

0.0025 

Table  1.  Bacterial  species  with  significantly  higher  abundance  in  Nod2~'~  HFD  mice  than  in  WT  HFD  mice*. 
*P<  0.05  for  all  species.  The  entire  microbiome  analysis  data  comparing  bacterial  diversity  in  the  stools  of  WT 
and  Nod2~'~  mice  maintained  on  HFD  and  LFD  have  been  deposited  in  NCBI  SRA,  accession  No.  SRP076031. 


of  pups  for  8  weeks.  Germ-free  pups  colonized  with  microbiota  from  Nod2~'~  HFD  mice  were  more  susceptible 
to  HFD-dependent  weight  gain  than  pups  colonized  with  microbiota  from  WT  HFD  mice,  as  manifested  by 
significantly  higher  weight  gain  (Fig.  7b).  Furthermore,  pups  colonized  with  Nod2~'~  HFD  microbiota  were  also 
more  susceptible  to  HFD-induced  metabolic  dysfunction,  as  they  had  significantly  higher  liver  cholesterol  and 
triglycerides,  and  higher  serum  glucose  than  pups  colonized  with  microbiota  from  WT  HFD  mice  (Fig.  7c-e), 
but  there  was  no  difference  in  serum  cholesterol  or  triglycerides  between  the  2  groups  of  mice.  In  conclusion, 
microbiota  from  Nod2~'~  HFD  mice  can  transfer  to  WT  outbred  germ-free  mice  enhanced  susceptibility  to 
HFD-dependent  weight  gain  and  metabolic  dysfunction. 

Discussion 

Multiple  factors,  including  host  genetics,  microbiota,  and  the  environment  regulate  sensitivity  to  obesity  and  met¬ 
abolic  disease  and  development  of  their  clinical  manifestations.  However,  the  specific  components  that  promote 
or  protect  from  these  diseases  are  poorly  understood.  There  is  increasing  evidence  for  cross  talk  between  innate 
immunity  and  metabolism  and  for  their  combined  role  in  the  development  of  obesity  and  associated  diseases6,26-28. 
However  there  is  no  previous  evidence  for  the  role  of  the  innate  immunity  gene  Nod2  in  regulating  sensitiv¬ 
ity  to  obesity.  Here  we  have  discovered  that  Nod2  and  Nod2- dependent  microbiota  protect  BALB/c  mice  from 
HFD-induced  pathology 

We  show  that  Nod2~/~  BALB/c  mice  on  HFD  become  obese  and  exhibit  many  key  manifestations  of  metabolic 
disease  including  hypercholesterolemia,  hypertriglyceridemia,  hyperglycemia,  and  hyperinsulinemia,  increased 
adiposity,  and  steatosis.  The  effect  of  Nod2  deletion  on  HFD-induced  obesity  is  strain-dependent,  as  Nod2~'~  and 
WT  C57BL/6  mice  on  HFD  gain  similar  weight,  but  Nod2~'~  C57BL/6  mice  exhibit  more  severe  insulin  resist¬ 
ance21  than  Nod2~'~  BALB/c  mice.  Thus,  Nod2  differentially  regulates  the  response  of  BALB/c  and  C57BL/6  mice 
to  HFD,  which  emphasizes  the  importance  of  the  hosts  entire  genetic  background  in  regulation  of  obesity  and 
metabolic  syndrome  and  reflects  the  variability  in  the  development  of  disease  in  humans. 

Metabolic  inflammation  plays  a  key  role  in  the  development  of  obesity  and  metabolic  diseases  in  humans 
and  in  animal  models  of  obesity5,6.  Our  results  also  demonstrate  that  HFD-dependent  obesity  in  Nod2~'~  mice 
correlates  with  increased  expression  of  a  large  number  of  genes  in  adipose  tissue  that  participate  in  immune 
responses.  They  include  genes  for  chemokines  and  cell  surface  molecules  involved  in  the  recruitment  and  activa¬ 
tion  of  T  cells,  monocytes,  dendritic  cells  and  eosinophils,  and  for  enzymes  in  eicosanoid  metabolism.  In  addi¬ 
tion,  Nod2~'~  HFD  mice  have  increased  infiltration  of  F4/80+  macrophages  and  enhanced  levels  of  the  transcripts 
of  Adgrel  (gene  for  F4/80)  and  CD68 ,  also  a  marker  for  monocytes/macrophages29,  in  the  adipose  tissue.  Similar 
to  our  results,  increased  inflammatory  markers  were  also  observed  in  Nod2~'~  C57BL/6  mice  on  HFD21,  which 
suggests  that  inflammation  likely  plays  an  important  role  in  the  development  of  metabolic  dysfunction  in  both 
strains  of  mice. 

Development  of  obesity  is  usually  accompanied  by  dysfunction  in  lipid  and/or  carbohydrate  metabolism, 
and  our  results  also  demonstrate  that  HFD-dependent  obesity  in  Nod2^'~  mice  correlates  with  changes  in  the 
expression  of  many  genes  involved  in  intermediary  metabolism,  primarily  lipid  metabolism.  Several  of  these 
differentially  expressed  genes  and  their  products  have  an  established  role  in  the  development  of  obesity  or  are 
associated  with  metabolic  diseases30-35.  Leptin  is  an  adipokine  that  regulates  appetite  and  energy  balance,  and  in 
obese  individuals  there  is  decreased  sensitivity  to  leptin,  which  may  result  in  increased  appetite36.  Nod2~'~  mice 
on  HFD  have  elevated  levels  of  serum  leptin,  enhanced  transcription  of  leptin  ( Lep )  in  the  adipose  tissue,  and 
decreased  expression  of  the  leptin  receptor  in  the  liver  ( Lepr ).  These  results  suggest  that  Nod2~'~  mice  on  HFD 
develop  resistance  to  leptin,  however  there  was  no  difference  in  the  amount  of  food  consumed  between  WT 
and  Nod2~'~  mice  on  HFD.  Genes  known  to  participate  in  multiple  other  biological  pathways,  including  DNA 
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Figure  6.  Antibiotic  treatment  prevents  the  development  of  obesity,  steatosis,  and  metabolic  dysfunction  in 
Nod2~'~  mice  on  HFD.  (a)  Nod2~'~  and  WT  mice  were  treated  with  Abx  to  deplete  their  intestinal  microbiota. 
Control  mice  were  not  treated  with  antibiotics.  Three  weeks  after  the  start  of  Abx  all  mice  were  placed  on  HFD 
and  monitored  for  weight  gain,  histopathology,  and  metabolites,  (b)  Percent  gain  in  weight  compared  to  week  0 
of  HFD.  (c,d)  After  14  weeks  of  HFD,  mice  were  analyzed  for  liver  histopathology  in  BODIPY  and  DAPI  stained 
sections  and  (d)  percent  frequency  distribution  of  LD  area  (pm2)  for  -4,000  LDs/group.  (e-h)  After  14  weeks 
of  HFD,  (e)  liver  triglycerides,  (f)  liver  cholesterol,  (g)  serum  cholesterol,  and  (h)  serum  triglycerides  were 
measured.  (b,d)  Results  are  means  ±  SEM,  (C)  representative  images,  and  (e-h)  individual  data  of  6-10  mice/ 
group.  *P<  0.05  and  **P<  0.001,  Nod2~l~  Abx  +  HFD  versus  Nod2~'~  HFD. 


SCIENTIFIC  REPORTS  |  7:  548  |  DOI:10.1038/s41598-017-00484-2 


12 


www.nature.com/scientificreports/ 


a 


Pregnant 

WeekO 

germ-free 

Pups  born 

Week  3 

mice 

Mothers  and  pups 

Pups  weaned 

-1  week 

4,  3  weeks 

i 

5  weeks 

T 

T  T  T 

t  t 

t  T 

Gavage 
and  co-house 

Gavage  mothers 

Gavage  pups 

T 


HFD 


** 


WT  Nod2_/" 
HFD  microflora 


Figure  7.  Intestinal  microbiota  from  Nod2~'~  mice  on  HFD  increases  susceptibility  of  germ-free  mice  to 
HFD-induced  obesity  and  metabolic  dysfunction,  (a)  One  week  before  delivery,  we  colonized  germ-free  WT 
Swiss  Webster  female  pregnant  mice,  kept  under  sterile  conditions,  with  microbiota  from  WT  or  Nod2~'~ 
mice  maintained  on  HFD.  Colonization  was  done  by  both  gavaging  with  fecal  microbiota  and  by  co-housing. 
Gavaging  was  done  into  the  stomach  with  fecal  homogenates  from  WT  or  Nod2~'~  mice  maintained  on  HFD 
(indicated  by  j).  At  the  same  time,  the  gavaged  mice  were  placed  on  HFD  and  also  co-housed  with  Nod2~'~  or 
WT  female  mice  that  had  been  maintained  on  HFD  for  ~20  weeks.  The  co-housed  Nod2~'~  HFD  and  WT 
HFD  mice  were  removed  when  the  pups  were  born  (week  0).  Mothers  and  pups  were  maintained  on  HFD  for 
the  length  of  the  experiment,  (b-e)  Pups  were  monitored  for  (b)  gain  in  weight  every  week  and  after  8  weeks 
assayed  for  (c)  liver  cholesterol,  (d)  liver  triglycerides  and  (e)  serum  glucose,  (b)  Results  are  means  ±  SEM  and 
(c-e)  individual  data  of  19-21  mice/group.  *P<  0.05  and  **P<  0.001,  Nod2~'~  Abx+  HFD  versus  Nod2~'~ 
HFD. 


structure,  gene  expression,  cell  cycle,  development,  neural  processes,  transport,  cytoskeleton,  and  extracellular 
matrix  were  also  differentially  expressed  between  Nod2~'~  HFD  mice  compared  with  WT  HFD  mice,  which  sug¬ 
gest  the  development  of  a  complex  pathology 

Obesity  is  often  associated  with  non-alcoholic  fatty  liver  disease,  and  Nod2~/~  mice  on  HFD  also  exhibit 
similar  manifestations,  including  enlarged  fatty  liver,  elevated  levels  of  cholesterol  and  triglycerides  in  the  liver, 
and  large  LDs  in  the  hepatocytes.  LDs  are  composed  of  a  lipid  core  surrounded  by  a  phospholipid  monolayer  and 
are  coated  with  proteins  that  regulate  LD  physiology,  including  enzymes  in  lipid  metabolism  and  members  of 
the  perilipin  family,  which  regulate  lipid  metabolism  and  trafficking37,38.  The  development  of  LDs  is  regulated  by 
multiple  factors,  including  host  genetics,  composition  of  the  intestinal  microbiota,  and  diet37-40.  Using  proteomics 
we  identified  many  proteins  in  the  LD  enriched  fraction  with  significantly  changed  abundance  in  Nod2~'~  HFD 
mice  compared  with  WT  HFD  mice.  All  identified  proteins  have  known  functions  in  lipid  metabolism  and/or 
vesicle  trafficking  and  several  of  these  proteins,  for  example  Abhd5,  Pnpla3,  and  COPI,  have  an  established  role  or 
are  associated  with  obesity  and  /or  steatosis37,38,41-45.  In  addition,  many  genes  involved  in  lipid  synthesis,  degra¬ 
dation,  and  storage  were  differentially  regulated  in  the  liver  of  Nod2~'~  mice  on  HFD  compared  with  WT  mice  on 
HFD.  These  results,  thus,  provide  clues  for  the  molecular  mechanism  responsible  for  the  development  of  steatosis 
in  Nod2~'~  HFD  mice.  Moreover,  the  altered  microbiota  in  Nod2~'~  HFD  mice  also  contributes  to  increased  lipid 
accumulation  in  hepatocytes  and  may  play  a  role  in  the  formation  of  large  LDs. 

Obese  mice  and  humans  often  show  increased  abundance  of  Firmicutes  and  Proteobacteria ,  but  decreased 
abundance  of  Bacteroidetes3, 46, 47.  Accordingly,  our  results  also  show  higher  abundance  of  six  Firmicutes  species 
and  one  Proteobacteria  species  in  obese  Nod2~'~  HFD  mice  compared  with  non-obese  WT  HFD  mice.  The  abun¬ 
dance  of  several  of  these  species  is  also  increased  in  obese  humans48-51.  Our  results  comparing  the  abundance 
of  Bacteroidetes  in  obese  and  non-obese  mice  were  mixed:  three  Bacteroidetes  species  were  increased  and  two 
were  decreased  in  Nod2~'~  HFD  mice  compared  with  WT  HFD  mice,  consistent  with  many  exceptions  to  the 
general  rule  of  decreased  Bacteroidetes  in  obesity,  e.g.,  increased  abundance  of  Bacteroides  was  found  in  some 
obese  and  obese-prone  populations52, 53.  Also,  consistent  with  the  association  of  inflammation  with  obesity,  the 
abundance  of  two  members  of  Deferribacteres  were  increased  in  obese  Nod2~'~  HFD  mice  compared  with  all 
other  non-obese  groups  of  mice,  as  Deferribacteres  are  early  markers  of  inflammation54.  Our  results  concur  with 
previous  data  showing  dysbiosis  in  Nod2~'~  C57BL/6  mice  on  HFD21,  but  the  specific  differences  were  not  the 
same  as  those  seen  in  Nod2~'~  BALB/c  mice  on  HFD. 
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There  is  increasing  evidence  indicating  a  role  for  the  gut  microbiota  in  the  development  of  obesity  in  humans 
and  in  animal  models  of  obesity2, 3.  Microbiota  with  increased  capacity  for  harvesting  energy,  referred  to  as  “obese 
microbiota”  has  been  identified  in  obese  individuals,  and  furthermore,  the  obese  microbiota  from  humans  can 
transfer  the  obesity  phenotype  to  mice55, 56 .  Our  results  showing  abrogation  of  obesity  and  metabolic  changes  in 
Nod2~'~  HFD  mice  by  antibiotics  and  transfer  of  increased  sensitivity  to  obesity  in  germ-free  mice  by  microbi¬ 
ota  from  Nod2~'~  HFD  mice  suggest  that  Nod2  may  protect  mice  from  harboring  such  “obese  microbiota”  This 
notion  is  further  supported  by  a  decreased  total  energy  content  in  the  stools  of  Nod2~/~  HFD  mice  compared 
with  WT  HFD  mice,  suggesting  greater  extraction  of  energy  from  food  in  Nod2~'~  HFD  mice,  leading  to  obesity 
and  metabolic  disease.  These  results  indicate  that  there  are  many  similarities  in  the  mechanisms  responsible  for 
microbiota-induced  obesity  in  humans  and  mice,  although  the  specific  bacteria  present  in  humans  and  mice  may 
differ. 

In  summary,  we  have  discovered  that  deletion  of  Nod2  abolishes  the  resistance  of  BALB/c  mice  to 
HFD-induced  obesity  and  metabolic  dysfunction.  The  development  of  obesity  and  metabolic  dysfunction  in 
Nod2~'~  mice  on  HFD  has  many  similarities  to  the  development  of  HFD-induced  pathology  in  humans,  which 
makes  Nod2~'~  BALB/c  mice  a  relevant  model  for  the  study  of  HFD-induced  obesity  and  associated  diseases. 

Methods 

Mice.  Nod2~'~  mice  on  BALB/c  background  were  described  previously  and  deletion  of  the  Nod2  genes  was 
confirmed  by  PCR  analysis  of  genomic  DNA19,20.  The  original  founder  WT  BALB/c  mice  were  obtained  from 
Harlan-Sprague-Dawly.  WT  germ-free  mice  (Swiss  Webster  female  and  male)  were  obtained  from  Taconic  Farms 
(Hudson,  NY).  All  WT  and  Nod2~'~  mice  on  BALB/c  background  were  bred  and  kept  under  conventional  path¬ 
ogen-free  conditions  in  the  same  room  in  our  facility  to  minimize  the  influence  of  differences  in  the  environ¬ 
ment.  For  each  experiment,  mice  from  several  different  cages  and  breeder  pairs  were  used.  To  avoid  changes  in 
microbiome  that  could  accumulate  over  extended  period  of  time,  we  backcross  our  mutant  mice  to  WT  females 
once  every  other  year  and  re-derive  our  homozygous  knockout  breeding  pairs.  The  latter  strategy  also  minimized 
genetic  drift  in  the  population.  Some  investigators  use  knockout  and  WT  littermates  from  heterozygous  breeding 
pairs  to  minimize  parent-to-parent  and  cage-to-cage  variations.  We  do  not  use  this  strategy  for  two  reasons:  first, 
this  strategy  may  skew  the  results  to  the  particular  microbiota  present  only  in  this  breeding  pair.  Second  and 
most  importantly,  we  noticed  that  the  effect  of  innate  immunity  genes  on  the  composition  of  microbiota  is  not 
instantaneous,  but  takes  time,  and  stabilization  of  microbiota  characteristic  of  a  given  knockout  may  take  more 
than  one  generation20.  The  BALB/c  background  of  knockout  mice  and  their  negative  status  for  all  common  viral 
and  bacterial  pathogens  and  parasites  were  confirmed  as  previously  described19.  The  Indiana  University  School 
of  Medicine-Northwest  Institutional  Animal  Care  and  Use  Committee  approved  all  experiments  with  mice.  All 
methods  were  performed  in  accordance  with  the  relevant  guidelines  and  regulations. 

Obesity  model  in  conventional  mice.  WT  and  Nod2~'~  female  mice  in  BALB/c  background,  matched 
for  weight  and  age  (6-7  weeks  old)  were  randomly  divided  into  two  groups  and  started  on  HFD  (fat  60kcal%, 
D12492)  or  LFD  (fat  10kcal%,  D12450B)  from  Research  Diets  and  maintained  on  the  special  diets  until  the  end 
of  the  experiment  (14  or  30  weeks).  The  total  calories/g  of  HFD  and  LFD  were  the  same.  Mice  were  weighed 
every  week  and  the  amount  of  food  consumed  each  week  was  measured  to  confirm  that  the  differences  in  obe¬ 
sity  are  not  due  to  differences  in  food  intake.  One  to  two  weeks  before  sacrifice,  stool  samples  were  collected  for 
microbiota  analysis,  bomb  calorimetry  and  for  microbiota  transplant  experiments,  fasting  glucose  was  assayed, 
and  glucose  tolerance  and  insulin  resistance  assays  performed.  At  the  time  of  sacrifice,  blood  was  collected  for 
biochemical  assays;  visceral  adipose  tissue  and  liver  were  weighed  and  subcutaneous  adipose  tissue  was  measured 
with  a  micrometer.  Visceral  adipose  tissue  and  liver  were  processed  for  histology,  RNA,  and  flow  cytometry.  Lipid 
droplets  were  purified  from  the  liver  and  proteins  identified  by  mass  spectrometry  We  focused  on  adipose  tissue 
and  liver,  because  they  are  the  major  organs  involved  in  intermediary  metabolism. 

Obesity  in  antibiotic-treated  mice.  We  depleted  the  microbiota  in  WT  and  Nod2~'~  mice  using  an  estab¬ 
lished  antibiotic  treatment20, 23“25.  Mice  were  matched  for  weight  and  age  and  randomly  divided  into  two  groups. 
One  group  of  4  weeks  old  mice  was  started  on  an  antibiotic  mix  (Abx)  of  ciprofloxacin  (lOOpg/ml)  and  metro¬ 
nidazole  (250  pg/ml)  in  their  drinking  water.  After  three  weeks  of  Abx  treatment  anaerobic  stool  cultures  were 
negative  (<  103  bacteria/g  of  feces,  freshly  collected  in  reducing  buffer),  and  mice  were  placed  on  HFD  and  Abx 
was  continued,  and  body  weight  and  food  consumed  were  measured  each  week.  Control  mice  were  started  on 
HFD  at  the  same  age,  but  not  treated  with  Abx.  At  14  weeks  after  the  start  of  HFD,  mice  were  sacrificed,  blood  was 
collected  for  biochemical  assays  and  liver  was  processed  for  histology  and  lipid  assays. 

Obesity  in  microbiota  transplanted  germ-free  mice.  One  week  before  delivery,  we  colonized 
germ-free  WT  Swiss  Webster  female  pregnant  mice,  kept  under  sterile  conditions,  with  microbiota  from  WT 
or  Nod2~'~  mice  maintained  on  HFD.  Colonization  was  done  by  both  gavaging  with  fecal  microbiota  and  by 
co-housing.  Gavaging  was  done  into  the  stomach  with  10  mg  of  fecal  homogenates  (containing  5  x  109  bacteria 
in  0.2  ml)  prepared  from  freshly  defecated  stool  samples,  collected  from  WT  or  Nod2~'~  mice  maintained  on 
HFD  for  29-30  weeks,  and  immediately  placed  in  reduced  anaerobic  medium,  homogenized,  and  stored  frozen  at 
—80  °C,  as  previously  described20,25.  At  the  same  time,  the  gavaged  mice  were  placed  on  HFD  and  also  co-housing 
with  Nod2~'~  or  WT  female  mice  that  had  been  maintained  on  HFD  for  20  weeks  was  started  and  continued  until 
one  day  after  delivery  The  mothers  were  gavaged  weekly  and  maintained  on  HFD  until  weaning  of  the  pups  at  3 
weeks  of  age.  Pups  were  gavaged  at  weaning  and  once  a  week  thereafter,  and  continued  on  HFD.  Body  weight  was 
measured  every  week  and  blood  was  collected  at  the  time  of  sacrifice.  We  used  this  strategy  to  ensure  successful 
colonization  with  microbiota,  since  we  did  not  know  whether  gavaging  or  co-housing  would  be  more  effective 
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in  transferring  obesity-promoting  microbiota.  We  colonized  pregnant  mice  to  ensure  that  the  pups  will  be  born 
colonized  with  a  specific  gut  microbiota.  The  presence  of  microbiota  at  birth  is  essential  for  the  normal  develop¬ 
ment  of  host  metabolism  and  immune  systems  and  alterations  in  the  gut  bacteria  at  birth  increases  susceptibility 
to  diet-induced  obesity57. 

Histopathology  of  liver  and  adipose  tissue.  Liver  and  adipose  tissue  were  fixed  in  10%  buffered  forma¬ 
lin,  embedded  in  paraffin,  sectioned,  and  stained  with  hematoxylin  and  eosin.  Because  lipid  is  lost  with  paraffin 
embedding,  we  also  analyzed  frozen  liver  sections  stained  with  BODIPY  493/503  (Molecular  Probes)  and  DAPI 
(Vector  Laboratories).  The  liver  sections  were  fixed  in  4%  paraformaldehyde/PBS  for  10  min  at  room  tempera¬ 
ture,  rinsed  with  PBS,  and  stained  with  1  pg/ml  of  BODIPY  493/503  in  PBS  for  5  min  at  room  temperature.  The 
sections  were  rinsed  with  PBS  and  then  covered  with  Vectashield  mounting  medium  containing  DAPI.  The  area 
for  lipid  droplets  and  adipocytes  was  measured  using  ImageJ  and  Adiposoft58,  respectively 

Biochemical  analyses.  Mice  were  fasted  for  5h  and  blood  was  collected.  Serum  cholesterol  and  triglycer¬ 
ides  were  quantified  using  fluorometric  kits  from  Bio  Vision  Inc.  Serum  insulin  and  leptin  were  analyzed  using 
ELISA  kits  from  Linco  Research  Inc.  Blood  glucose  was  measured  using  the  CVS/pharmacy  Advanced  Glucose 
Meter  and  test  strips.  Cholesterol  and  triglycerides  were  also  quantified  in  liver  homogenates  (BioVision  Inc). 

Glucose  tolerance  test  and  estimation  of  insulin  sensitivity.  WT  and  Nod2~/~  mice  maintained  on 
HFD  for  29  weeks  were  fasted  for  5  h  and  baseline  glucose  measured  using  blood  collected  from  the  tail  vein.  For 
the  glucose  tolerance  test,  mice  were  injected  intrap eritoneally  with  2  g  glucose/kg  body  weight  in  sterile  PBS  and 
blood  glucose  measured  at  40,  80,  and  120  min  after  injection.  For  the  insulin  sensitivity  test,  mice  were  injected 
intraperitoneally  with  0.4U/kg  body  weight  of  Humulin-R  insulin  (Health Warehouse)  in  sterile  PBS  and  blood 
glucose  measured  at  15,  30,  60,  90,  and  120  min  after  injection. 

Bomb  calorimetry.  To  identify  changes  in  energy  extracted  from  the  food  between  Nod2~'~  and  WT  mice 
on  HFD,  fresh  fecal  samples  were  collected  from  individual  mice,  weighed,  and  frozen  at  —80  °C  until  processing. 
Samples  were  dried  at  58  °C  for  24  h,  dried  pellets  were  ground  and  energy  content  was  measured  using  a  Parr 
Oxygen  Bomb  Calorimeter  and  total  caloric  value  computed  as  per  the  manufacturers  instructions. 

Lipid  droplet  enrichment  and  protein  analysis.  Lipid  droplet  enrichment.  We  enriched  LDs  using  a 
modified  procedure59.  Freshly  dissected  livers  were  homogenized  on  ice  with  a  glass  Dounce  homogenizer  con¬ 
taining  ice-cold  buffer  A  (20  mM  tricine  and  250  mM  sucrose,  pH  7.8)  with  0.2  mM  PMSF  and  protease  inhibitor 
cocktail.  The  homogenate  was  centrifuged  at  100  g  for  10  min  at  4°C.  The  supernatant  was  further  disrupted  for 
15  min  at  500  psi  in  a  nitrogen  bomb  chamber  and  centrifuged  at  3,000  g  for  10  min  at  4  °C  to  remove  nuclei,  cell 
debris,  and  unbroken  cells.  The  post-nuclear  supernatant  (PNS)  was  transferred  into  a  15  ml  tube  and  3  ml  of 
buffer  B  (20  mM  HEPES,  100  mM  KC1  and  2  mM  MgCl2,  pH  7.4)  was  loaded  on  top  and  centrifuged  at  2,000  g 
for  30  min  at  4°C.  LDs  were  collected  from  the  top  band  and  washed  several  times  with  200  pi  buffer  B  and  cen¬ 
trifuged  at  20,000  g  for  5  min  at  4°C.  The  samples  were  mixed  with  chloroform/acetone  (0.25:0.75,  v/v),  kept  at 
— 80  °C  overnight,  and  proteins  and  lipids  were  separated  by  centrifugation  at  20,000  g  for  10  min  at  4°C.  The 
protein  pellet  was  air-dried,  dissolved  in  8  M  Urea  and  analyzed  by  mass  spectrometry 

Proteomic  analysis  by  mass  spectrometry  Proteins  in  the  LD  enriched  fraction  were  identified  by  mass  spectrom¬ 
etry  at  the  Laboratory  for  Biological  Mass  Spectrometry,  Indiana  University  Bloomington.  Proteins  were  digested 
with  trypsin  and  purified  peptides  were  run  on  an  Eksigent  Nano  2D  HPLC  coupled  to  a  LTQ  Velos  Pro  mass 
spectrometer  (Thermo  Fisher  Scientific,  Waltham  MA).  The  resulting  data  were  searched  in  Protein  Prospector 
against  the  Homo  sapiens  or  Mus  musculus  proteome  in  the  SwissProt  database.  Spectral  counting  was  used  for 
protein  quantification.  Normalized  data  for  comparison  were  obtained  by  dividing  spectral  count  of  a  single 
protein  by  the  sum  of  spectral  counts  in  the  same  sample.  The  fold  change  was  computed  for  Nod2~'~  HFD  with 
respect  to  WT  HFD  and  genes  with  a  fold  change  of  >1.5  or  <0.6  and  P<  0.05  were  considered  as  significantly 
different  between  the  2  groups  of  mice. 

RNA  transcriptomics.  For  identification  of  genes  that  were  differentially  expressed  between  Nod2~'~  and 
WT  mice  on  HFD,  we  analyzed  the  total  RNA  population  in  the  liver  and  adipose  tissue  by  RNA  sequencing. 
RNA  was  isolated  from  the  liver  and  adipose  tissue  of  individual  WT  and  Nod2~'~  mice  maintained  on  HFD 
for  30  weeks  using  the  TRIZOL  method  (InVitrogen),  followed  by  purification  on  RNeasy  spin  columns  using 
Qiagen  RNeasy  Mini  Kit60.  RNA  sequencing  and  initial  analysis  were  performed  at  the  Center  for  Genomics  and 
Bioinformatics,  Indiana  University  Bloomington.  98%  of  all  reads  (-311  million  reads)  passed  a  quality  filter  of 
q20  requiring  a  minimum  length  of  50  bases  after  adapter  and  quality  trimming,  94.7%  of  these  mapped  to  the 
mouse  genome  assembly  GRCm38.75  using  tophat  ver.2.0.10  and  95.5%  of  these  mapped  uniquely  The  reads 
were  annotated  using  the  software  DAVID61  and  differential  expression  with  statistical  analysis  was  calculated 
using  the  software  DESeq262.  The  fold  change  was  computed  for  Nod2~'~  HFD  with  respect  to  WT  HFD  and 
genes  that  had  a  fold  change  of  >2  or  <0.5  and  P  <  0.05  and  with  a  5%  FDR  were  considered  as  significantly 
different  between  the  2  groups. 

Isolation  of  stromal/vascular  cells  from  adipose  tissue  and  flow  cytometry.  Stromal/vascular 
fraction  from  adipose  tissue,  which  contains  immune  cells,  was  isolated  using  a  modified  procedure63.  Visceral 
adipose  tissue  was  isolated,  weighed,  and  digested  with  1  mg/ml  collagenase  D  and  0.15  pg/ml  DNAse  I  in  RPMI- 
5FCS  (RPMI-1640  supplemented  with  5%  FCS)  for  40  min  in  a  shaking  37  °C  water  bath  and  samples  were  cen¬ 
trifuged  at  1,000  rpm  for  10  min  at  room  temperature;  at  this  speed  adipocytes  stay  on  the  top  and  the  stromal/ 
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vascular  (SV)  cells  form  a  pellet.  The  supernatant  with  the  adipocytes  was  carefully  removed  and  discarded  and 
the  pellet  was  resuspended  in  RPMI-5FCS  and  filtered  through  a  100  [im  filter  to  remove  any  remaining  adi¬ 
pocytes.  SV  cells  were  stained  with  CD45-Vioblue  and  F4/80-APC  (Biolegend)  for  20  min,  analyzed  by  flow 
cytometry  using  MACSQuant  (Miltenyi)  cytometer  and  F4/80  positive  cells  were  measured  within  the  CD45  gate. 

Gut  microbiome  analysis.  Stool  collection  and  DNA  extraction.  We  collected  freshly  defecated  stools  from 
individual  female  WT  and  Nod2~'~  mice  (200  mg/mouse,  6  mice/strain)  maintained  for  28-30  weeks  on  HFD  or 
LFD  and  isolated  DNA  using  the  MoBio  Fecal  DNA  isolation  kit25.  For  each  genotype,  mice  originated  from  three 
different  litters  from  different  parents  and  2  mice  per  litter),  weaned  into  separate  cages,  and  all  cages  were  kept  in 
the  same  room  in  our  animal  facility  This  strategy  allows  stabilization  of  microbiota  and  minimizes  the  variability 
observed  between  different  litters  due  to  different  parents  and  different  cages  and  rooms20,25,60. 

Pyrosequencing  of  16S  rRNA  genes  and  microbiome  analysis.  Genetic  phylotyping  (community  profiling) 
using  pyrosequencing  of  the  variable  regions  of  bacterial  16  S  ribosomal  RNA  (rRNA)  genes  with  16  S  univer¬ 
sal  Eubacterial  primers  530  F  (5'-GTG  CCA  GCMGCN  GCG  G)  and  1100R  (5'-GGG  TTN  CGN  TCG  TTG) 
was  performed  using  Roche  454  technology  and  bTEFAP®  pipeline  at  MR  DNA  Molecular  Research  LP25, 64.  A 
single-step  30  cycle  PCR  using  HotStarTaq  Plus  Master  Mix  Kit  (Qiagen)  was  used  under  the  following  condi¬ 
tions:  94  °C  for  3  min,  followed  by  28  cycles  of  94  °C  for  30  sec,  53  °C  for  40  sec,  and  72  °C  for  1  min,  and  elonga¬ 
tion  step  at  72  °C  for  5  min.  Following  PCR,  all  amplicon  products  from  different  samples  were  mixed  in  equal 
concentrations  and  purified  using  Agencourt  Ampure  beads  (Agencourt  Bioscience  Corporation).  Samples 
were  sequenced  using  Roche  454  FLX  titanium  instrument  and  reagents  by  following  manufacturer  s  guidelines. 
Sequencing  data  were  processed  using  MR  DNA  Molecular  Research  LP  analysis  pipeline  (www.mrdnalab.com), 
including  depletion  of  barcodes  and  primers,  removal  of  short  sequences  (<200  bp),  sequences  with  ambigu¬ 
ous  base  calls,  and  sequences  with  homopolymer  runs  exceeding  6  bp,  de-noising,  and  removal  of  chimeras65. 
Operational  taxonomic  units  (OTUs)  were  then  defined  after  removal  of  singleton  sequences  and  clustering, 
and  taxonomically  classified  using  BLASTn  against  a  curated  GreenGenes  database  and  assigned  to  species  at 
>97%  identity  to  the  reference  sequence  (<3%  divergence),  to  genus  at  95%-97%  identity,  family  at  90%-95% 
identity,  order  at  85%-90%  identity,  class  at  80%-85%  identity,  phylum  at  77%-80%  identity,  and  unclassified  at 
<77%  identity66.  We  compared  the  sequences  and  analyzed  the  changes  in  bacterial  ecology  using  Quantitative 
Insights  Into  Microbial  Ecology  (QIIME;  http://qiime.sourceforge.net),  Bayesian,  PyNAST,  and  UniFrac67-69.  The 
Shannon  diversity  index  H  =  -Xp2  In (p2)  and  Shannon  equitability  index  EH  =  H/ln(S),  where  p{  is  the  proportion 
of  the  ith  taxonomic  unit  (OTU  or  species)  and  S  is  the  total  number  of  taxonomic  units,  were  calculated  using 
Microsoft  Excel.  The  average  total  number  of  reads  per  mouse  was  6978,  with  no  significant  differences  between 
the  numbers  of  reads/mouse  between  WT  and  Nod2~'~  and  HFD  and  LFD  mice,  and  detected  bacteria  with  a 
frequency  of  >0.001%. 

Statistical  analysis.  The  significance  of  differences  between  the  numbers  of  recovered  bacterial  species  and 
OTUs  was  determined  by  the  z-test  using  Microsoft  Excel.  The  significance  of  differences  between  all  other  quan¬ 
titative  results  presented  as  means  ±  SEM  were  determined  by  the  Students  f-test  using  Microsoft  Excel  or  by 
ANOVA  using  Prism.  The  N  and  P  values  are  indicated  in  the  Figures  and  Tables;  P  <  0.05  was  considered  signif¬ 
icant.  The  heatmaps  were  generated  using  Java  TreeView  and  represent  mean  fold  changes  in  Nod2~'~  HFD  mice 
relative  to  WT  HFD  mice,  after  converting  <1  ratios  to  negative  fold  difference  using  the  formula:  (—1)/  ratio. 
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Supplementary  Figure  SI.  Nodi1'  and  WT  mice  on  HFD  show  no  difference  in  the  amount  of  food  consumed 
and  Nodi1'  and  WT  mice  on  regular  chow  show  no  difference  in  voluntary  activity,  (a)  The  amount  of  food 
consumed  by  WT  and  Nodi1'  mice  on  HFD  per  cage  was  measured  each  week  for  30  weeks  and  the  amount  of 
food  consumed  per  mouse  is  shown.  Results  are  means  of  4  cages  for  each  group,  and  total  number  of  mice, 
A/=20-24.  (b)  Nodi1'  and  WT  mice  on  regular  chow  have  no  difference  in  voluntary  wheel  activity. 


MA  plots  for  adipose  genes  in  Nod2~^~  and  WT  HFD  mice 
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MA  plots  for  liver  genes  in  Nod2~^~  and  WT  HFD  mice 
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Supplementary  Figure  S2.  Nodi /_  mice  on  HFD  have  many  differentially  expressed  genes  in  the  liver  and 
adipose  tissue  compared  with  WT  mice  on  HFD.  MA  plots  showing  gene  expression  in  the  (A)  liver  and  (B) 
adipose  tissue  of  Nodi1'  HFD  mice  relative  to  WT  HFD  mice.  Each  dot  represents  a  gene.  The  Y-axis  is  the  log2 
fold  change  ( Nodi1 '  HFD  /WT  HFD)  and  the  X-axis  is  the  average  of  the  counts  normalized  by  size  factor.  All 
genes  falling  on  the  red  horizontal  line  intersecting  at  0  on  the  Y-axis  have  a  log2  ratio  of  zero  (no  change  in 
gene  expression).  All  the  genes  falling  above  0  are  up  regulated  in  Nodi1' ,  while  those  below  0  are  down 
regulated.  The  genes  colored  red  are  significantly  up  or  down  regulated  (significance  defined  at  P< 0.05  at  5% 
FDR).  The  triangles  represent  genes  that  are  differentially  expressed  over  log2  fold  change. 
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Supplementary  Figure  S3,  a-diversity  in  stool  microbiota  in  WT  HFD,  NodZ1'  HFD,  WT  LFD,  and  NodZ /_  LFD 
mice,  (a)  Total  numbers  of  species,  numbers  of  species/mouse,  Shannon  diversity  index,  and  Shannon 
equitability  index  for  bacterial  species,  (b)  Total  numbers  of  OTUs,  numbers  of  OTUs/mouse,  Shannon  diversity 
index,  and  Shannon  equitability  index  for  bacterial  OTUs.  The  results  are  means  ±  SEM  or  totals;  N  =  6 
mice/group;  *P< 0.05  WT  versus  NodZ1'  (on  HFD  or  LFD);  #P< 0.05  HFD  versus  LFD  (for  WT  or  NodZ1').  The 
entire  microbiome  analysis  data  comparing  bacterial  diversity  in  the  stools  of  WT  and  NodZ1'  mice  maintained 
on  HFD  and  LFD  have  been  deposited  in  NCBI  SRA,  accession  No.  SRP076031. 
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*  Clostridia  ■  Candidatus  saccharibacteria  W  Alphaproteobacteria 

■■  Bacteroidia  ■  Deltaproteobacteria  ■  Cyanobacteria 

^  Gammaproteobacteria  ■  Flavobacteriia  w  Spirochaetia 

r  Bacilli  ■  Actinobacteria  Deferribacteres 

r  Epsilonproteobacteria  ■  Erysipelotrichia  ^  Betaproteobacteria 
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Supplementary  Figure  S4.  Class  (J-diversity  in  stool  microbiota  in  WT  HFD,  NodZ1'  HFD,  WT  LFD,  and  NodZ1' 
LFD  mice.  Class  abundance  expressed  as  %  of  total  stool  microbiota  is  shown  for  individual  mice  (1-6)  and  an 
average  (Avg)  for  all  bacterial  classes  (top  graphs)  and  for  17  less  abundant  classes  expressed  as  %  of  total  low 
abundant  classes  (excluding  two  most  abundant  classes,  Clostridia  and  Bacteroida)  (bottom  graphs);  *,  classes 
with  significantly  (at  P< 0.05)  increased  abundance  in  WT  versus  NodZ1'  mice  (on  HFD  or  LFD);  #,  classes  with 
significantly  (at  P< 0.05)  increased  abundance  in  HFD  versus  LFD  groups  (for  WT  or  NodZ1'  mice).  The  entire 
microbiome  analysis  data  comparing  bacterial  diversity  in  the  stools  of  WT  and  NodZ1'  mice  maintained  on 
HFD  and  LFD  have  been  deposited  in  NCBI  SRA,  accession  No.  SRP076031. 


